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STABILITY OF VISCOUS ST. VENANT ROLL-WAVES: FROM ONSET TO 

INFINITE-FROUDE NUMBER LIMIT 


BLAKE BARKER, MATHEW A. JOHNSON, PASCAL NOBLE, L.MIGUEL RODRIGUES, AND KEVIN ZUMBRUN 

Abstract. We study the spectral stability of roll-wave solutions of the viscous St. Venant equations 
modeling inclined shallow-water flow, both at onset in the small-Froude number or “weakly unstable” 
limit F —> 2 + and for general values of the Froude number F, including the limit F —> + 00 . In the former, 

F —> 2 + , limit, the shallow water equations are formally approximated by a Korteweg de Vries/Kuramoto- 
Sivashinsky (KdV-KS) equation that is a singular perturbation of the standard Korteweg de Vries (KdV) 
equation modeling horizontal shallow water flow. Our main analytical result is to rigorously validate this 
formal limit, showing that stability as F —¥ 2 + is equivalent to stability of the corresponding KdV-KS 
waves in the KdV limit. Together with recent results obtained for KdV-KS by Johnson-Noble-Rodrigues- 
Zumbrun and Barker, this gives not only the first rigorous verification of stability for any single viscous St. 

Venant roll wave, but a complete classification of stability in the weakly unstable limit. In the remainder 
of the paper, we investigate numerically and analytically the evolution of the stability diagram as Froude 
number increases to infinity. Notably, we find transition at around F = 2.3 from weakly unstable to 
different, large-F behavior, with stability determined by simple power law relations. The latter stability 
criteria are potentially useful in hydraulic engineering applications, for which typically 2.5 < F < 6.0. 


1. Introduction 

In this paper, we investigate the stability of periodic wavetrain, or roll-wave, solutions of the inclined 
viscous shallow-water equations of St. Venant, appearing in nondimensional Eulerian form as 

(1.1) d t h + d x (hu ) = 0, d t (hu) + d x ^ hu 2 + = h - \u\u + vd x {hd x u), 

where F is a Froude number, given by the ratio between (a chosen reference) speed of the fluid and speed 
of gravity waves, and u = R~ x , with R e the Reynolds number of the fluid. System (1.1) describes the 
motion of a thin layer of fluid flowing down an inclined plane, with h denoting fluid height, u fluid velocity 
averaged with respect to height, x longitudinal distance along the plane, and t time. The terms h and 
|it| u on the righthand side of the second equation model, respectively, gravitational force and turbulent 
friction along the bottom. 1 

Roll-waves are well-known hydrodynamic instabilities of (1.1), arising in the region F > 2 for which 
constant solutions, corresponding to parallel flow, are unstable. They appear in the modeling of such 
diverse phenomena as landslides, river and spillway flow, and the topography of sand dunes and sea beds; 
see Fig. 1 (a)-(b) for physical examples of roll waves and Fig. 1 (c) for a typical wavetrain solution of (1.1). 
As motivated by these applications, their stability properties have been studied formally, numerically, 
and experimentally in various physically interesting regimes; see, for example, [BM04] for a useful survey 
of this literature. However, up until now, there has been no complete rigorous stability analysis of viscous 
St. Venant roll-waves either at the linear (spectral) or nonlinear level. 

Recently, the authors, in various combinations, have developed a theoretical framework for the study of 
nonlinear stability of these and related periodic waves. Specifically, for the model at hand, it was shown 
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in [JZN11] that, under standard diffusive spectral stability assumptions (conditions (D1)-(D3) in §1.1.1) 
together with a technical “slope condition” ((1.4) below) satisfied for “moderate” values 2 < F fz 3.5 of 
F and a genericity assumption ((HI) below) satisfied almost everywhere in parameter space, 2 roll-waves 
are nonlinearly stable in the sense that localized perturbations converge to localized spatial modulations 
of the background periodic wave. See also [BJN + 12, BJN + 13] for discussions in the related context of 
the Kuramoto—Sivashinsky equation [Kur84, KT75, Siv77, Siv83]. More recently, for general (partially) 
parabolic systems, detailed nonlinear asymptotic behavior under localized and nonlocalized perturbations 
has been established in [JNRZ14] in terms of certain formal modulation, or “Whitham,” equations. 3 

This reduces the study of stability and asymptotic behavior, at least for moderate values of F. to 
verification of the spectral stability conditions (D1)-(D3), concerning Floquet spectrum of the associated 
eigenvalue ODE. However, it is in general a hard problem to verify such spectral assumptions analytically. 
Indeed, up to now, spectral stability has not been rigorously verified for any roll wave solution of the 
viscous St. Venant equations (1.1). 



Figure 1. Roll waves (a) on a spillway and (b) in the lab: pictures courtesy of Neil 
Balmforth, UBC. (c) Periodic profile of (1.1), F = \/6, u = 0.1, q = 1.5745, X = 17.15. 

For better comparison to experiment, we extended the profile here as constant in transverse 
direction. 

In some particular situations, for example, at the onset of hydrodynamical instability, analytical proof of 
spectral stability may be possible using perturbation techniques. However, most of the known examples 
concern reaction diffusion equations and related models like the Swift Hohenberg equations, Rayleigh 
Benard convection or Taylor Couette flows that are all described, near the instability threshold of a 
background constant solution, by a Ginzburg-Landau equation derived as an amplitude equation [Mie02], 
Associated with classical Hopf bifurcation, this normal form may be rigorously validated in terms of 
existence and stability by Lyapunov-Schmidt reduction about a limiting constant-coefficient operator 
[CE90, Mie97a, Mie97b]. 

By contrast, the corresponding model for onset of hydrodynamic (roll-wave) instability in (1.1) is, at 
least formally, the Korteweg-de Vries/Kuramoto-Sivashinsky equation (KdV-KS) 

(1.2) dsv + vdyv + edyv + 5(dyv + dyv) = 0, V5 > 0, VF G 1, 

with 0 < i « 1, £ > 0, a singular perturbation of the Korteweg-de Vries (KdV) equation 4 . Equation 
(1.2) is derived as an amplitude equation for the shallow water system (1.1) near the critical value F = 2 
above which steady constant-height flows are unstable, in the small-amplitude limit h = h + 5 2 v and in 
the KdV time and space scaling (V, S) = (6(x — cot),5 3 t) with 6 = \JF — 2, where Co is an appropriate 
reference wave speed: see Section 2.1 below for details in the Lagrangian formulation. Alternatively, it 
may be derived from the full Navier-Stokes equations with free boundary from which (1.1) is derived in 

Indeed, this appears numerically to be satisfied for all profiles. 

3 See [OZIO] for the easier multidimensional case, in which behavior is asymptotically linear due to faster decay of the 
linearized propagator. 

4 Without loss of generality, one can assume that e 2 + 8 2 = 1. See [BJN + 13]. 
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the shallow-water limit, for Reynolds number R near the critical value R c above which steady Nusselt 
flows are unstable; see [Win93, YY03]. 

In this case, neither existence nor stability reduce to computations involving constant-coefficient op¬ 
erators; rather, the reference states are arbitrary-amplitude periodic solutions of KdV, and the relevant 
operators (variable-coefficient) linearizations thereof. This makes behavior considerably richer, and both 
analysis and validation of the amplitude equations considerably more complicated than in the Ginzburg- 
Landau case mentioned above. Likewise, onset occurs for (1.1) not through Hopf bifurcation from a single 
equilibrium, but through Bogdanev-Takens, or saddle-node bifurcation involving collision of two equilib¬ 
ria, as discussed, e.g., in [HC87, BJRZ11], with limiting period thus +oo, consistent with the 1/5 spatial 
scaling of the formal model. (The standard unfolding of a Bogdanev-Takens bifurcation as a perturbed 
Hamiltonian system is also consistent with KdV-KS; see Remark 2.1.) 

Nevertheless, similarly as in previous works by Mielke [Mie97a, Mie97b] in the reaction diffusion setting, 
where the stability of periodic waves for the amplitude (Ginzburg-Landau) equation provides a stability 
result for periodic waves of the full system (Swift Hohenberg equation or Rayleight Benard convection), 
we may expect that stability for the amplitude equation, here the KdV-KS equation (1.2) will provide 
some information on the stability of periodic waves for the viscous St. Venant system (1.1), at least 
in the weakly-unstable limit F —> 2 + . Our first main goal is to rigorously validate this conjecture, 
showing that stability of roll-waves in the weakly unstable limit F —> 2 + is determined by stability of 
corresponding solutions of (1.2) under the rescaling described above. Together with previous results 
[BN95, JNRZ15, Barl4] on stability for (1.2), this gives the first complete nonlinear stability results for 
roll waves of (1.1): more, it gives a complete classification of stability in the weakly unstable limit. 

This gives at the same time a rigorous justification in a particular instance of the much more generally 
applicable and better-studied (1.2) as a canonical model for weak hydrodynamic instability in inclined 
thin-ffim flow; see, e.g., [BN95, CD02, CDK93, PSU07]. Looked at from this opposite point of view, (1.1) 
gives an interesting extension in a specific case of (1.2) into the large-amplitude, strongly-unstable regime. 
Our second main goal is, by a combination of rigorous analysis and (nonrigorous but numerically well- 
conditioned) numerical experiment, to continue our analysis into this large-amplitude regime, performing 
a systematic stability analysis for F on the entire range of existence F > 2 of periodic roll-wave solutions 
of (1.1). Our main finding here is a remarkably simple power-law description of curves bounding the 
region of stability in parameter space from above and below, across which particular high-frequency 
and low-frequency stability transitions occur. These curves eventually meet, yielding instability for F 
sufficiently large. The large-T description is quite different from the small-F description of weakly 
unstable theory; indeed, there is a dramatic transition from small- to large-F behavior at F ~ 2.3, with 
behavior governed thereafter by the large-F version. This distinction appears important for hydraulic 
engineering applications, where F is typically 2.5 — 6.0 and sometimes 10 — 20 or higher [Jef25, Bro69, 
Bro70, AeM91, RG12, RG13, FSMA03]. Far from the onset, nonlinear stability for spectrally stable waves 
remains an interesting open question due to failure of the technical (slope) condition (1.4); as discussed 
below, there is reason to think this may be dropped. 

1.1. Summary of previous work. We begin by recalling some known results that will be relied upon 
throughout our analysis. In particular, we begin by recalling how spectral stability (in a suitable diffusive 
sense) may provide a detailed nonlinear stability result, a fact that strongly underpins and motivates our 
spectral studies. We then recall the relevant numerical and analytical results for the amplitude equation 
(1.2), upon which our entire weakly-unstable analysis for 0<L-2<1 hinges. 

1.1.1. Diffusive spectral stability conditions. We first recall the standard diffusive spectral stability condi¬ 
tions as defined in various contexts in, for example, [Sch98, Sch96, JZ11, JZ10, JZN11, BJN + 13, JNRZ14]. 

Let u(x, t ) = u(x — ct ) define a spatially periodic traveling-wave solution of a general partial differential 
equation dtu = F{u) with period (without loss of generality) one, or, equivalently, a be a stationary 
solution of dtu = F + cd x u with period one, and let L := (dF / du){u) + cd x denote the associated 
linearized operator about u. As L is a linear differential operator with 1-periodic coefficients, standard 
results from Floquet theory dictate that non-trivial solutions of Lv = Xv can not be integrable on M, 
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more generally they can not have finite norm in L P (R) for any 1 < p < oo. Indeed, it follows by standard 
arguments that the T 2 (M)-spectrum of L is purely continuous and that A £ a L 2 m\(L) if and only if the 
spectral problem Lv = \v has an L°° (M)-eigenfunction of the form 

v(x; A,£) = e* x w( x;A,£) 


for some £ £ [tt, tt) and w £ L 2 ei .([0,1]); see [Gar93] or [Rodl3, p.30-31] for details. In particular, 
A £ a L 2 ^(L) if and only if there exists a £ £ [—7r,7r) such that there is a non-trivial 1-periodic solution 
of the equation 


L^w = A w, where 


(L^w) (x) := e~* x L 



(x). 


and 

c7 l 2 (ir)(-C') = = [J a L2 er ([ 0jl ])(^)- 

£S[—7T,7r) 


The parameter £ is referred to as the Bloch or Floquet frequency and the operators L^ are the Bloch 
operators associated to L. Since the Bloch operators have compactly embedded domains in L 2 er ([0, T]) 
their spectrum consists entirely of discrete eigenvalues that depend continuously on the Bloch parameter 
£. Thus, the spectrum of L consists entirely of L°° (R)-eigenvalues and may be decomposed into countably 
many curves A(£) such that A(£) £ cr L2 er ([o,i])(-k£) for £ £ [— tt, tt). 

Suppose, further, that u is a transversal 5 orbit of the traveling-wave ODE F(u) + cd x u = 0. Then 
near u, the Implicit Function Theorem guarantees a smooth manifold of nearby 1-periodic traveling-wave 
solutions of (possibly) different speeds, with some dimension N £ N, 6 not accounting for invariance under 
translations. Then, the diffusive spectral stability conditions are : 


(Dl) <t l2(r) (L) C{A | SKA < 0} U {0}. 

(D2) There exists a 0 > 0 such that for all £ £ [— tt, tt) we have ir([0J]) (L c )c{A|»A<-^| 2 }. 
(D3) A = 0 is an eigenvalue of Lq with generalized eigenspace £o C L 2 er ([0, 1]) of dimension N. 

Under mild additional technical hypotheses to do with regularity of the coefficients of IF, hyperbolic- 
parabolic structure, etc., conditions (D1)-(D3) have been shown in all of the above-mentioned settings 
in particular for periodic waves of either (1.1) or (1.2) — to imply nonlinear modulation stability, at 
Gaussian rate : more precisely, provided || (?i — zZ) |i=o 11 Z/i (R)ni^ s (M) is sufficiently small for some s sufficiently 
large, there exists a function ip(x,t) with ^(x,0) = 0 such that the solution satisfies 


(1.3) Ju(qt)-u(--V’M)-ct)|| £ p ( R) + ||V*,t^(-,i)||LP(R) < C(l + t) 2 (1 1 /p\ 2 < p < oo, 

valid for all t > 0; see [JZ11, JZ10, JZN11, BJN + 13, JNRZ14]. In the case of (1.1), (1.2), for which 
coefficients depend analytically on the solution, essentially there suffices the single technical hypothesis: 
(HI) The N zero eigenvalues of Lq split linearly as £ is varied with |£| sufficiently small, in the sense 
that they may be expanded as A_,-(£) = oy£ + o(£) for some constants ay £ C distinct. 

We note that, since the existence of the expansion Ay(£) = ay£ + o(£) in (HI) may be proved indepen¬ 
dently, the hypothesis (HI) really concerns distinctness of the ay, which is equivalent to the condition that 
the characteristics of a (formally) related first-order Whitham modulation system be distinct, a condition 
that, in the case of analytic dependence of the underlying equations, as here, either holds generically with 
respect to nondegenerate parametrizations of the manifold of periodic traveling waves, or else uniformly 
fails. For (1.1), there is an additional slope condition 

(1.4) h x /h < (cuFY 2 , 


used to obtain hyperbolic-parabolic damping and high-frequency resolvent estimates by Kawashima-type 
energy estimates, necessary to obtain the desired nonlinear modulational stability result; see [JZN11, 
Section 4.3]. Condition (1.4) is known to be sufficient but a priori not necessary. As discussed in Section 


J In a sense compatible with the algebraic structure of the system. 

5 For both (1.1) and (1.2), an easy dimensional count gives N = 2 [JZN11, BJN + 10] because of the presence of one local 
conservation law in the respective sets of equations. 
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1.3, we believe that this can be replaced by its average, 0 < {cvF)~ 2 : hence dropped 7 . However condition 
(1.4) holds evidently in the small-amplitude limit F —> 2 + , and is observed numerically for moderate 
values 2 < F < 3.5. 

The above nonlinear stability results motivate a detailed analytical inspection of the conditions (Dl)- 
(D3) and (HI), which is precisely the intent of the weakly nonlinear analysis presented in Section 2 
below. We note that these conditions may be readily checked numerically, in a well-conditioned way, using 
either Hill’s method (Galerkin approximation), or numerical Evans function analysis (shooting/continuous 
orthogonalization); see [BJN + 12, BJN + 13, BJN + 10]. 

1.1.2. Numerical evaluation for viscous St. Venant and KdV-KS. The diffusive spectral stability condi¬ 
tions (D1)-(D3) have been studied numerically for the viscous St. Venant equations (1.1) in [BJN + 10] 
for certain “typical” waves and Froude numbers F , with results indicating existence of both stable and 
unstable waves: more precisely, the existence of a single “band” of stable waves as period is varied for 
fixed F. This echoes the much earlier numerical study of roll-waves of the classical Kuramoto-Sivashinsky 
equation (KS) (e = 0 for (1.2)) in [FST 86 ] and elsewhere, that obtained similar results. 

Equations (1.2) have received substantially more attention, as canonical models for hydrodynamical 
instability in a variety of thin-film settings; as derived formally in [CDK93], see also [Rodl3, p.16, 
footnote 10], the model (1.2) with the addition of a further term D(vy)y, D constant, gives a general 
form for such instabilities in the weakly unstable regime. A systematic numerical study of this more 
general model was carried out in [CDK93], across all values of e, 6, D, and the period X of the wave, and, 
by different methods in [BJN + 13], for the value D = 0 only; see Figure 2, reprinted from [BJN + 13] (in 
close agreement also with the results of [CDK93]). As noted in [CDK93], it may be observed from Figure 
2 that the small stable band for e/5 <C 1 enlarges with addition of dispersion/decrease in 5, reaching its 
largest size at S/e = 0 (corresponding to the singular KdV) limit. For intermediate ratios of 5/e, behavior 
can be considerably more complicated, with bifurcation to multiple stable bands as this ratio is varied. 

KdV-KS stability boundaries 
301 
25 
20 
x 15 
10 
51 


0 0.2 0.4 0.6 0.8 

e 

Figure 2. Stability boundaries (in period X ) vs. parameter e for the KdV-KS equation 
( 1 . 2 ) with e 2 + 5 2 = 1. 



1.1.3. The KdV limit 5 0 + . Of special interest for us is the KdV limit 5 -A 0 + for (1.2), treated with 

varying degrees of rigor in [EMR93, BN95, NR13, JNRZ14, Barl4]: a singularly perturbed Hamiltonian 
indeed, completely integrable — system. We cite briefly the relevant results; for details, see [JNRZ14, 
Bar 14]. 

Proposition 1.1 (Existence [EMR93]). Given any positive integer r > 1, there exists So > 0 such that 
there exist periodic traveling wave solutions v$(6), 6 = Y — a$S, of (1.2) (with e = 1) that are analytic 


7 See the Authors Note at the end of Section 1.3 for recent progress in this direction. 
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functions of 9 £ M and C r functions of 5 £ [0,c>o). When r > 3, profiles vs expand as 6 
2-parameter family 

( v s (9; a 0 , k ) = T 0 (9; a 0 , k, k) + 5Ti(9) + 5 2 T 2 (9) + 0(5 3 ), 

\ a s = a 0 (a 0l k,K) + 5 2 cr 2 +0(5 3 ), 


0 + 


as a 


(1.5) 
where 

To(9;ao,k,K) = ao + 12k 2 K 2 cn 2 (k 9, k ), <to = ao + 4 k 2 (2k 2 — 1), 
comprise the 3-parameter family (up to translation) of periodic KdV profiles and their speeds; cn(-,&:) is 
the Jacobi elliptic cosine function with elliptic modulus k £ (0,1); ao is a parameter related to Galilean 
invariance; and k = G(k) is determined via the selection principle 


(K(k)Q(k)\ 


7 


V 


7r 


) 


2(k A - k 2 + 1 )E(k) - (1 - k 2 )(2 - k 2 )I<(k) 


20 (-2 + 3k’ 2 + 3A; 4 - 2k 6 )E(k) + (k 6 + k 4 - 4 k 2 + 2 )K(k) ’ 

where K(k ) and T(/c) are the complete elliptic integrals of the first and second kind. The period X(k ) = 
2K(k)/Q(k) is in one-to-one correspondence with k. Moreover the functions (T))j = i )2 are (respectively 
odd and even) solutions of the linear equations 


( 1 . 6 ) 


£ 0 [T 0 ]Ti = T" + T"", £ 0 [T 0 ]T 2 =[ J j--a 2 T 0 ) + T" + T"", 


on (0,2 I\(k)/Q(k)) with periodic boundary conditions, where ToPo] := — — dg (Tq — ao) denotes the 

linearized KdV operator about To- 

Throughout this manuscript, we let ^(-; ao, k) and < 7 , 5 ( 00 , k) denote the periodic traveling wave profiles 
and wave speeds, respectively, as described in Proposition 1.5. 


Remark 1.2. The additional parameter n for periodic KdV waves as compared to KdV-KS waves reflects 
the existence of the additional conserved quantity of the Hamiltonian at 5 = 0. The selection principle 
k = Q(k) is precisely the condition that the periodic Hamiltonian orbits at 5 = 0 persist, to first order, 
for 0 < 5 <C 1. Alternatively, this condition can be written explicitly as j 2K W)/ K 7 ^ (Tq + T()")dx = 0. 

By Galilean invariance of the underling equation (1.2), the stability properties of the above-described 
X = X(fc)-periodic solutions are independent of the parameter ao- Hence, for stability purposes one may 
identify waves with a common period, fixing a 0 and studying stability of a one-parameter family in k. It 
is known [KSF84, Spe 88 , BD09] that the spectra of the linearized operator -Cpo], considered on L 2 (M), 
about a periodic KdV wave To is spectrally stable in the neutral, Hamiltonian, sense, i.e., all eigenvalues 
of the Bloch operators 

£ ? [T 0 ] := (dr + if) {-(dy + if) 2 -T 0 + a 0 ) : L 2 pei ( 0, X) -> L 2 er (0, X), 

considered with compactly embedded domain TTp er (0, X ), are purely imaginary for each f £ [— 7 t/X, 7 t/X). 
Moreover, the explicit description of the spectrum obtained in [BD09] also yields 8 that A = 0 is an 
eigenvalue of -CoPo] of algebraic multiplicity three, that A = 0 is an eigenvalue of -Ccpo] only if f = 0, 
and that the three zero eigenvalues of £oPo(s k, G(k))\ expand about for |£| < 1 as 

(1.7) A j(f) = iatj(f)f = ifaj + 0(f 2 ), j = 1,2,3 with oij £ M distinct. 

We introduce one final technical condition, first observed then proved numerically to hold, at least for 
KdV waves that are limits as <5 —> 0 of stable waves of (1.2) [BD09, BJN + 13, Barl4 ]: 9 

(A) A given parameter k £ (0,1) is said to satisfy condition (A) if the non-zero eigenvalues of the lin¬ 
earized (Bloch) KdV operator -C^po] about To(-; ao, k,Q(k)) are simple for each f £ [— ir/X, n/X). 

8 In [BJN + 13, Barl4], some of these facts remained unnoticed to the authors. In particular, in [BJN + 13], the condition 
that only £ = 0 yields A = 0 was gathered to condition (A) below to form condition (Al) and distinctness of the cq was 
identified as condition (A2). 

9 In fact (A) has been verified (see Proposition 1.4 below) on essentially the entire range k £ (0,1); we know of no instance 
where it fails. 
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Note that the set of k G (0,1) for which property (A) holds is open. 

Given a periodic traveling wave solution To(--,ao,k,Q(k)) of the KdV equation with elliptic modulus 
k G (0,1) satisfying condition (A) above, we now consider the spectral stability of the associated family 
of periodic traveling wave solutions vg(-',ao,k), defined for 5 G [0, <5o), with 5q as in Proposition 1.5, as 
solutions of the KdV-KS equation (1.2). To this end, notice that, assuming k G (0,1) satisfies condition 
(A), the non-zero Bloch-eigenvalues A(£) of the linearized KdV-KS operator 

Hv s ] = e~* [-<5 (4 + dy) ~ dy ~ dy (v s - a s )] e* : L 2 peT (0, X) L 2 er (0, X) 

admit a smooth 10 expansion in S for 0 < <5 <C 1. In particular, for each pair (£, Ao) with Ao G cr(£^[To])\{0} 
and £ G [— n/X, it/X) there is a unique spectral curve A(£, Ao ,S) bifurcating from Ao smoothly in <5, and 
it takes the form 

(1-8) A(<5; £, A 0 ) = A 0 + <5Ai(£, A 0 ) + 0(6 2 ) 

for some Ai(£, Ao). It is then natural to expect that the signs of the real parts of the first order correctors 
Ai(£, A 0 ) in the above expansion be indicative of stability or instability of the near-KdV profiles us for 
0 < <5 <C 1. With this motivation in mind, for any k G (0,1) that satisfies condition (A) above 11 , we 
define 

(1.9) Ind(A’) := sup 3? (Ai(£, Ao)). 

A o e<7(££[To(-;ao,fc,0(fc))])\{O} 

te[-*/X{k),n/X(k)) 

Evidently, Ind(fc) > 0 is a sufficient condition for the spectral instability for 0 < 5 <C 1 of the near-KdV 
waves vs bifurcating from To. The next proposition states that the condition Ind(fc) < 0 is also sufficient 
for the diffusive spectral stability of the vs, for 0 < 5 <C 1. Define the open set 

V := { k G (0,1) | condition (A) holds for k and Ind(fc) < 0 } . 

Proposition 1.3 (Limiting stability conditions [JNRZ15]). For each k G V there exists a neighborhood 
Dfc C (0,1) of k and 5o(k) > 0 such that fC V and for any ( ao,k,S ) G R x x (0,e>o(A;)) the non¬ 
degeneracy and spectral stability conditions (HI) and (D1)-(D3) hold for vs(-',do,k). In particular, V is 
open and do(-) can be chosen uniformly on compact subsets ofV. 

Proposition 1.3 reduces the question of diffusive spectral stability and nonlinear stability — in the sense 
defined in Section 1.1.1 — of the near-KdV profiles constructed in Proposition 1.5 to the verification of 
the structural condition (A) and the evaluation of the function Ind(fc). Note condition (A) is concerned 
only with the spectrum of the linearized KdV operator about the limiting KdV profile no; its validity is 
discussed in detail in [JNRZ14, Section 1], Further, note that due to the triple eigenvalue of the KdV 
linearized operator at the origin, the fact that Ind(fc) is sufficient for stability is far from a foregone 
conclusion, and this represents the main contribution of [JNRZ15]. 

To evaluate Ind(fc), using the complete integrability of the KdV equation to find an explicit parametriza- 
tion of the KdV spectrum and eigenfunctions about vq [BD09], one can construct a continuous multi¬ 
valued mapping 12 

[—7r /X, tt/A) xRi3 (£, Ao) 5R(Ai(£, A 0 )) G R 

that is explicitly computable in terms of Jacobi elliptic functions; see [BN95] or [BJN + 13, Appendix A.l]. 
This mapping may then be analyzed numerically. Numerical investigations of [BN95, BJN + 13] suggest 
stability for limiting periods X = X{k ) in an open interval (X m ,XM) and instability for X outside 
[A m , Am], with X m ps 8.45 and Xm ~ 26.1, thus completely classifying stability for 0 < 8 <C 1. The 
following result of Barker [Bar 14], established by numerical proof using interval arithmetic, gives rigorous 

10 In the sense of Proposition 1.1 that one can reach arbitrary prescribed regularity. 

^Condition (A) is independent of ao, holding for every ao or for none. Likewise, Ind(fc) is independent of ao- 
12 While Ai(£, Ao) is defined above only when Ao is a simple eigenvalue of the limiting linearized KdV operator Lj[uo] , it 
possesses an explicit expression in terms of k and an auxiliary Lax spectral parameter that extend this function to points 
where simplicity fails. Note in particular that this extension is triple-valued at (0,0). 
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validation of these observations for all limiting periods X(k ) except for a set near the boundaries of the 
domain of existence X(k) G (27r, +oo) ~ (6.2832, +oo) corresponding to the limits k —> 0,1. 

Proposition 1.4 (Numerical proof [Barl4]). With k m j n = 0.199910210210210 and k max = 0.999999999997, 
corresponding to X m m ~ 6.284 and X max ~ 48.3, condition (A) holds on [fc m i n , A: max ]. Moreover, there 
are ki G [0.9421,0.9426] and k r G [0.99999838520,0.99999838527], corresponding to Xi G [8.43,8.45] 
and X r G [26.0573,26.0575], such that V n [fc m i n , fe max ] = ( ki,k r ) and Ind takes positive values on 

[Gnin • k max ] \ [ki,k r \. 

As noted in [Barl4], the limits k —> 0 and k —>• 1 not treated in Proposition 1.4, corresponding to Hopf 
and homoclinic limits, though inaccessible by the numerical methods of [Bar 14], should be treatable by 
asymptotics relating spectra to those of (unstable [BJRZ11, BJN + 10]) limiting constant and homoclinic 
profiles. 


1.2. Description of main results. As mentioned previously, the fact that the KdV-KS equation (1.2) 
serves as an amplitude equation for the shallow water system (1.1) in the weakly nonlinear regime 
0 < 6 = \JF — 2 <C 1 suggests that Proposition 1.3 and Proposition 1.4 should have natural counterparts 
for the stability of roll-waves in (1.1). To explore this connection, following [JZN11, JNRZ14, BJRZ11], 
we find it convenient to rewrite the viscous St. Venant equations (1.1) in their equivalent Lagrangian 
form: 

(1.10) d t T - d x u = 0, d t u + d x 1 = 1-tu 2 + vd x {r~ 2 d x u), 

where r := 1/h and x denotes now a Lagrangian marker rather than a physical location x, satisfying the 
relations dt/dx = u(x,t) and dx/dx = r(x,t). In these coordinates, slope condition (1.4) takes the form 

(1.11) 2vu x < F~ 2 . 

Hereafter, we will work exclusively with the formulation (1.10). 


Remark 1.5. Though nontrivial, the one-to-one correspondence between periodic waves of the Eulerian 
and Lagrangian forms is a well-known fact. A fact that seems to have remained unnoticed until very 
recently is that this correspondence extends to the spectral level even in its Floquet-by-Floquet descrip¬ 
tion. In particular without loss of generality one may safely study spectral stability in either formulation. 
See for instance [BGNR14] for an explicit description of the former correspondence and [BGMR] for the 
spectral connection, both in the closely related context of the Euler-Korteweg system. 


1.2.1. The weakly unstable limit F —> 2 + . Our first three results, and the main analytical results of this 
paper, comprise a rigorous validation of KdV-KS as a description of roll-wave behavior in the weakly 
unstable limit F -> 2 + . Let (tq,uq), uq = 7o~ 1//2 , be a constant solution of (1.10), and cq := r^/2. 
Setting S = \[F — 2, we introduce the rescaled dependent and independent variables 


( 1 . 12 ) 


TO 


f = 3e> 2 - 1 , u = 65 2 i - 1 , Y = 


u 0 


_ Tq /4 5(x - c 0 t) 


,1/2 


S := 


5 3 t 


4tq/ 4 iA/ 2 


Our first result concerns existence of small amplitude periodic traveling wave solutions in the limit 6 —> 0 + . 


Theorem 1.6 (Existence). There exists 5o > 0 such that there exist periodic traveling wave solutions 
of (1.10), in the rescaled coordinates (1.12) (t,u)s{0), 6 = Y — cs(ao,k)S, that are analytic functions of 
0 G R, (do, k) eRx (0,1) and 5 G [0, do) and that, in the limit 5 —»• 0 + , expand as a 2-parameter family 

j" f s (9 ; a 0 , k) = ~Vg(0; a 0 , k ) + 0(5 2 ), 

us (O', d 0 , k) = —fs(9; d 0 , k) + 4f(<?(d 0 , k) - 3r 0 _1 C5(d 0 , k)f s (6; a 0 , k)), 

c<s(d 0 , k) = crj(d 0 , k) + 0(5 2 ), 
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(1.13) 


< 



where 5 = 5/(2 t^ V 1 / 2 ), v$(0;ao,k ) and a$(ao,k) are the small-5 traveling-wave profiles and speeds of 
KdV-KS described in (1.5) and 

q(a 0 ,k) = 24fe 2 (l - k 2 )(G(k)) 4 - a 0 (|a 0 + 4(Q (k)) 2 (2k 2 - 1)) = v\f2 -gqVq + v'q 

is a constant of integration in the limiting KdV traveling-wave ODE (see (2.7)). 

For brevity, throughout the paper we shall often leave implicit the dependence on (ao, k) or a o- 

Remark 1.7. As we will see in the analysis below, the weakly unstable limit for the St. Venant equations 
(1.10) is a regular perturbation of KdV, rather than a singular perturbation as in the KdV-KS case, a 
fact reflected in the stronger regularity conclusions of Theorem 1.6 as compared to Proposition 1.1. 

Our next result concerns the spectral stability of the small amplitude periodic traveling wave solutions 
constructed in Theorem 1.6 when subject to arbitrary small localized (i.e. integrable) perturbations on 
the line. 

Theorem 1.8 (Limiting stability conditions). For each k G V there exists a neighborhood D& C (0,1) of 
k and 5o(k) > 0 such that for any (ao, k, 5) G M x D*, x (0, 5o(k)) the non-degeneracy and spectral stability 
conditions (HI) and (D1)-(D3) hold for (r,u)s(-;ao,k), along with slope condition (1-11). In particular, 
5o(-) can be chosen uniformly on compact subsets ofV. Conversely for each k G (0,1) such that condition 
(A) holds but Ind(k) > 0 (where Ind is defined as in (1.9),), there exists a neighborhood Ilk V (0,1) of k 
and 5 0 (k) > 0 such that if (ao,k,6) G R x Ilf. x (0,5o(k)) then (r, u)g(--, ao, k) is spectrally unstable. 

From Theorem 1.6, it follows in particular that our roll waves have asymptotic period ~ and 
amplitude ~ 5 2 in the weakly unstable limit F —> 2 + ; that is, this is a long-wave, small-amplitude 
limit. Likewise, u x r-j 5 3 so that (1.11) is automatically satisfied for 5 <C 1. In Theorems 1.6 and 1.8, 
rescaling period and amplitude to order one, we find that in this regime existence and stability are indeed 
well-described by KdV-KS —> KdV: to zeroth order by KdV, and to first correction by KdV-KS. 

Combining Theorem 1.8 with Proposition 1.4, and untangling coordinate changes, we obtain the fol¬ 
lowing essentially complete description of stability of viscous St. Venant roll-waves in the limit F —> 2 + . 

Corollary 1.9 (Limiting stability region). For 5 = \/F — 2 sufficiently small, uniformly for 5X on 

compact sets, periodic traveling waves of (1.10) are stable for (Lagrangian) periods X G u 5/4 (A;,A r ) 

V 5 

1/2 1/2 

and unstable for X G -[A min , A/ and X G J W-(V r , V max ] where A min , X h X r , A max are as in 

V 5 v s 

Proposition l.f, in particular, X mw ~ 6.284, Xi ps 8.44, X r ps 26.1, and A max ps 48.3. 

1.2.2. Large-Froude number limit F — > +oo. We complement the above weakly nonlinear analysis by 
continuing into the large amplitude regime, beginning with a study of the distinguished large-Froude 
number limit F —> +oo. The description of this limit requires a choice of scaling in the parameters 
indexing the family of waves. To this end, we first emphasize that a suitable parametrization, available 
for the full range of Froude numbers, is given by (q, X), where q := —cf — u is a constant of integration 
in the associated traveling-wave ODE, corresponding to total outflow , and X is the period. As discussed 
in Section 3.1 below, the associated two-parameter family of possible scalings may be reduced by the 
requirements that (i) the limiting system be nontrivial, and (ii) the limit be a regular perturbation, to a 
one-parameter family indexed by a > — 2, given explicitly via 

(1.14) r = aF a , u = bF ~ a / 2 , c = CoF" 1 - 3 "/ 2 , A = X 0 F~ l / 2 ~ 5a / A , q = q 0 F~ a / 2 . 

where a, b : M —>• R and Co, Ao,® are re& l constants. Note, from the relation A = 1/k between period 
and wave number k, that we have also k = k$F 1 / 2 + 5a / 4 . 

Under this rescaling, moving to the co-moving frame (x,t) i—)• ( k(x — ct),t), we find that A'-periodic 
traveling wave solutions of (1.10) correspond to Ao-periodic solutions of the rescaled profile equation 

(1.15) a" = (—a 2 /®^^) (^oo / F’ _3 / 2_3a / 4 (co — 1/a 3 ) — 1 + a(® — ® F~ l a) 2 — 2cokQv(a') 2 / a?), 

where b is recovered from a via the identity b = —® — coF~ l a] see Section 3.1 below for details. Noting 
that the behavior of p-3/2~3a/A a s F —>• oo depends on whether a = —2 or a > —2, one finds two 
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different limiting profile equations in the limit F — > oo: a (disguised 13 ) Hamiltonian equation supporting 
a selection principle, when a > —2, and a non-Hamiltonian equation in the boundary case a = —2; 
see Section 3.1 below. Further, by elementary phase plane analysis when a > —2 or direct numerical 
investigation when a = — 2, periodic solutions of the limiting profile equations are seen to exist as 2- 
parameter families parametrized by the period Xq and the discharge rage q o- Noting that, by design, the 
rescaled profile equation (1.15) is a regular perturbation of the appropriate limiting profile equation for 
all a > —2, we readily obtain the following asymptotic description 

Theorem 1.10. For sufficiently large F, generically, X^-periodic profiles of (1.15), obtained under the 
scaling (1.14), emerge for each a > — 2 from X^-periodic solutions of the appropriate limiting profile 
equation obtained by taking F —» oo in (1.15) and, when a = —2, satisfying a suitable selection principle. 

Next, we study the spectral stability of a pair of fixed periodic profiles (a, b) constructed above. One 
may readily check that, under the further rescaling Fb = b and ^ ll / 2 + Q / 4 \ = A, the linearized spectral 
problems around such a periodic profile (a, b) is given by 

A a — cokoa' — k^b' = 0 

F~ 3/ 2 -3a/4 — cokob' — ko(a/d 3 )'^j = —2 F~ l abb — b 2 a + vkl(b'a 2 + 2coa a/a 3 )', 

where (a, 6) denotes the perturbation of the underlying state (a, b). The limiting spectral problems ob¬ 
tained by taking F —> oo again depend on whether a = —2 or a > —2. In particular, we note the spectral 
problem is Hamiltonian, and hence possesses a natural four-fold symmetry about the real- and imaginary- 
axes, when a > —2; see 3.1 below for details. Since (1.16) is, again by design, a regular perturbation of the 
appropriate limiting spectral problems for all a > —2, we obtain by standard perturbation methods (e.g., 
the spectral/Evans function convergence results of [PZ04]) the following sufficient instability condition. 

Corollary 1.11. For all a > —2, under the rescaling (1.14), the profiles of (1.15) converging as F —> oo 
to solutions of the appropriate limiting profile equation, as described in Theorem 1.10 are spectrally 
unstable if the appropriate limiting spectral problem about the limiting profiles admit L 2 (M) -spectrum in 
A with positive real part. 

As clearly discussed in Section 3.1, for a > —2 the limiting profile equation, associated selection 
principle, and limiting spectral problem are independent of the value of a. Thus, the above instability 
criterion for F —> +oo can be determined by the study of just two model equations: one for a = —2 
and one for any other fixed a > —2. Both regimes include particularly physically interesting choices 
since a = 0 corresponds to holding the outflow q constant as F —> oo, while a = —2 corresponds to 
holding the Eulerian period E(Xo) constant as F —> oo 14 . In Section 3.2 we investigate numerically the 
stability of the limiting spectral problems in both the cases a = —2 and a = 0. This numerical study 
indicates that, in both these cases, all periodic solutions of the appropriate limiting profile equations are 
spectrally unstable and hence, by Corollary 1.11, that spectrally stable periodic traveling wave solutions 
of the viscous St. Venant system (1.10) do not exist for sufficiently large Froude numbers; see Figure 6. 

Numerical Observation 1. For both a > —2 and a = —2, the limiting eigenvalue equations have 
strictly unstable spectra, hence converging profiles are spectrally unstable for F sufficiently large. 

1.2.3. Intermediate F. We complete our stability investigation in Section 3.3 by carrying out a numerical 
study for F bounded away from the distinguished values 2 and +oo of the L 2 (M)-spectrum of the linearized 
operator obtained from linearizing (1.10) about a given periodic traveling wave solution. For F relatively 
small (2 < F < 4), we find, unsurprisingly, a smooth continuation of the picture for F —> 2 + , featuring 
a single band of stable periods between two concave upward curves; see Fig. 4 (c). However, continuing 
into the large-but-not-infinite regime (2.5 < F < 100), we find considerable additional structure beyond 
that described in Numerical Observation 1. 

1 o 1 

Indeed, the profile equation in this case is Hamiltonian in the unknown h = - 

14 Here, the Eulerian scaling relation E = Eo F~^~^ follows by 5 = J ( ' Y f(x)dx, (1.14), and convergence of the profile a. 
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Figure 3. In (a) and (b) we plot a numerical sampling of the (unstable) spectrum cor¬ 
responding to the F —> oo limiting spectral problems for the cases a = —2 and a > —2, 
respectively, for a representative periodic stationary solution of the appropriate limiting 
profile equation. 


Namely, for a E [—2,0] we see that the stability region is enclosed in a lens-shaped region between 
two smooth concave upward curves corresponding to the lower stability boundary and an upper high- 
frequency instability boundary, pinching off at a special value F*(a) after which, consistent with Numerical 
Observation 1, stable roll waves no longer exist; see Figure 4(a) for the case a = —2. Examining these 
curves further for different values of a, we find that they obey a remarkably simple power-law description 
in terms of F. q , a, and X; see Figure 4(b) for an example log-log plot in the case a = —2. The general 
description of this power-law behavior is provided by the following. 

Numerical Observation 2. Both lower and upper stability boundaries appear for F 1 to be 
governed by universal power laws c\ log F+c^ log q+cs log * = d, independent of parameters —2 < a < 0, 
u > 0, where for the lower boundary, c\ = 0.692, C2 = —3.46, C3 = 1, and d = 0.3, and for the upper 
boundary, c\ = 0.791, C2 = —1.73, C3 = 1, and d = 6.22: see Figure 7. Values a > 0 were not computed. 

Together with the small-F description of Theorem 1.8, these observations give an essentially complete 
description of stability of periodic roll wave solutions of the St. Venant equations (1.1), for —2 < a < 0. 
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Figure 4. Lower and upper stability boundaries for a = —2, v = 0.1, and, motivated by 
(1.14), scaling q = 0AF~ 2 . Solid dots show numerically observed boundaries. Pale dashes 
indicate approximating curves given by (a) (upper) X = e 0 ' 08 ' E 1 ’ 88 and (lower) X = 
e - 2- 97 _p2-83 ) (b) (upper) log(X) = 1.881og(F) + 0.087 and (lower) log(X) = 2.831og(F) — 
2.97. Pale dotted curves (Green in color plates) indicate theoretical boundaries as F —> 2 + . 
(c) Small- to large-F transition. 


1.3. Discussion and open problems. There have been a number of numerical and analytical studies 
of viscous roll waves in certain small-amplitude limits, in particular for the KdV-KS equations governing 
formally the weakly unstable limit F —> 2 + [FST86, CD02, CDK93, BN95, EMR93, PSU07, BJN + 13, 
JNRZ15, Barl4], However, to our knowledge, the present study represents the first systematic investi¬ 
gation of the stability of arbitrary amplitude roll wave solutions of the viscous St. Venant equations for 
inclined thin film flow. 
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Figure 5. Convergence to Dressier waves: We plot f(x) with pale solid curves and f'(x) 
with dark dashed curves (green and blue respectively in color plates). Here a = —2, 
v = 0.1, q = 0AF~ 2 , and (a) F = 5, X & 6.25, (b) F = 10, X « 33.3, (c) F = 15, 
X « 107, (d) F = 5, X « 20.8, (e) F = 10, X « 83.3, (f) F = 15, X w 205. 


Our main mathematical contribution is the rigorous validation of the formal KdV-KS —> KdV limit 
as a description of behavior in the small Froude number/weakly unstable limit F —> 2 + . This, together 
with the works [EMR93, BN95, JNRZ15, Barl4] on KdV-KS—>-KdV, gives a complete classification of 
existence and stability of viscous St. Venant roll waves in the weakly unstable regime. We note again that 
KS-KdV—>KdV is a canonical weakly unstable limit for the type of long-wave instabilities arising in thin 
film flow, in the same way that the (real and complex) Ginzburg-Landau equations are canonical models 
for finite-wavelength “Turing-type” instabilities. However, its analysis, based on singular perturbations 
of periodic KdV solutions, is essentially different from that of the hnite-wavelength case based on regular 
perturbation of constant solutions [Mie97a, Mie97b, Sch96]. 

From a practical point of view, the main point is perhaps the numerically-obtained description of 
behavior in the passage from small-amplitude to large-amplitude behavior. In particular, the universal 
scaling law of Numerical Observation 2 gives an unexpected global, simple-to-apply description of stability 
that seems potentially of use in biological and engineering applications, for which the St. Venant equations 
appear to be the preferred ones in current use. (Compare with the very complicated behavior in Fig. 2 
as 5 is varied away from the small-d limit.) This adds new insight beyond the qualitative picture afforded 
by the canonical KdV-KS—>KdV limit. In particular, our numerical results indicate a sharp transition 
at F ~ 2.3 from the quantitative predictions of the small-amplitude theory to the quite different large-F 
prediction of Numerical Observation 2. As hydraulic engineering applications typically involve values 
2.5 < F < 20 [Jef25, Bro69, Bro70, AeM91, RG12, RG13, FSMA03], this distinction appears physically 
quite relevant. 

A very interesting open problem, both from the mathematical and engineering point of view, is to 
rigorously verify this numerically-observed rule of thumb. As noted above (see Fig. 4), the upper and lower 
stability boundaries described in Numerical Observation 2 obey different scalings from those prescribed 
in (1.14), as F —> oo, with period growing faster than X ~ ^?-i/2-5«/4 ^y a f ac t or i? 1 / 2 -ci+a(c2/2+5/4) 
that is 1 for a < ct* ~ 1.54: in particular, for the two main physical values of interest a = —2 and 
a = 0, corresponding to constant (Eulerian) period and constant inflow, respectively. Indeed, given the 
large values of F to which the stability region extends, this may be deduced by Numerical Observation 1, 
which implies that all such waves of period 0(F -1 / 2 ~ 5 °7 4 ) are necessarily unstable for F>1. 

An important consequence is that, rescaling viscosity v so that the resulting period Xv after standard 
invariant scaling remains constant, we find that u —> 0. Hence, the limiting behavior of stable waves is 
described by the joint inviscid, large-Froude number limit v —> 0, F —> +oo. Figure 5, depicting periodic 
profiles at the upper and lower stability boundaries for values F = 5,10,15, clearly indicate convergence 
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as F increases to inviscirl Dressier waves [Dre49], alternating smooth portions and shock discontinuities. 
This agrees with recent observations of [BL05] that large-amplitude roll waves are experimentally well- 
predicted by a simplified, asymptotic version of the inviscid theory. We conjecture that our lower (low- 
frequency) stability boundary, corresponding to loss of hyperbolicity of associated Whitham equations, 
agrees with the inviscid threshold suggested by [Nob06, Theorem 1.2], 15 while the upper (high-frequency) 
stability boundary, corresponding to appearance of unstable spectra far from the origin, arises through a 
homoclinic, or “large-V,” limit similar to that studied in [Gar97, SS01] for reaction diffusion, KdV and 
related equations. 

We note that both analysis and numerics are complicated in the large-V limit by the appearance, 
differently from the case treated in [SS01] of essential spectra through the origin of the limiting solitary 
wave profile at X = +oo, along with the usual zero eigenvalue imposed by translational invariance. 
Whereas point spectra of a solitary wave are approximated as X —»• +oo by individual loops of Floquet 
spectra, curves of continuous spectra are “tiled” by arcs of length ~ V -1 , leading to the plethora of 
zero-eigenvalues (marked as pale dots, red in color plates) visible in Fig. 8(c)—(d). The large number 
of roots as X —> oo leads to numerical difficulty for both Hill’s method and numerical Evans function 
techniques, making the resolution of the stability region an extremely delicate computation, requiring 40 
days on IU’s 370-node Quarry supercomputer cluster to complete (Appendix C). The asymptotic analysis 
of this region is thus of considerable practical as well as theoretical interest. 

Another very interesting open question in the large-F regime is nonlinear stability of spectrally stable 
wave, in the absence of slope condition (1.4), or 2 vu x < F~ 2 , where u denotes the background profile. 
This condition was used in two ways in the nonlinear stability analysis of [JZN11]: first, to carry out 
high-frequency resolvent estimates used to obtain linearized stability estimates, and second, to carry out 
“nonlinear damping” type estimates used to close the nonlinear iteration yielding nonlinear modulational 
stability. In [JZN11], both were carried out by closely-related energy estimates. However, in the linear 
case, applying more delicate linear ODE techniques like those of Section 2.2, one finds that the needed 
high-frequency resolvent bounds require not that (1.11) hold pointwise but only on average. As u is 
periodic, the resulting averaged condition is 0 < F~ 2 , hence always satisfied. Thus, condition (1.11) 
seems superfluous for linearized stability. This strongly suggests that (1.11) may be dropped in the 
nonlinear analysis as well, either by a pseudodifferential analysis paralleling that of the linear case, or by 
an energy estimate with strategically-chosen (periodic) exponential weight analogously as in [Zum07] for 
viscous shocks. 

Finally, recall that, just as the KdV-KS—>KdV limit is derived formally from the viscous St. Venant 
equations in the weakly unstable limit, the viscous St. Venant equations are derived formally from 
the more fundamental free-boundary Navier-Stokes equations in the shallow water limit. Alternatively, 
the KdV-KS—7-KdV limit may be derived directly from the Navier-Stokes equations in a formal weakly 
unstable/shallow water limit. Rigorous verification of this formal limit, directly from the free-boundary 
Navier-Stokes equations, is perhaps the fundamental open problem in the theory. 

1.4. Plan of the paper. In Section 2.1, we recall the formal derivation of KdV-KS from the viscous St. 
Venant system and establish Theorem 1.6 concerning the existence of small amplitude roll-waves. These 
calculations will serve as a guideline for the subsequent analysis: indeed, we will follow the general strategy 
for the proof of spectral stability for KdV-KS periodic waves presented in [JNRZ15]. In Section 2.2, we 
begin studying the stability of these small amplitude roll-waves by computing a priori estimates on possible 
unstable eigenvalues for their associated linearized (Bloch) operators: energy estimates provide natural 
0(1) bounds (as <5 = \JF — 2 —> 0) whereas an approximate diagonalization process is needed to obtain 
the sharper bound 0(<5 3 ). At this stage, we recover, after a suitable rescaling and up to some negligible 
terms, the spectral problem associated to KdV-KS obtained after the Fenichel’s transformations. In 
Section 2.3, we then follow the proof in [JNRZ15] to complete the spectral stability analysis: for any 
fixed non-zero Bloch number, possible unstable eigenvalues A(£) for the linearized St. Venant system as 
A — >• 0 + are expanded as A(<5; £, Ao) = <5 3 Ao + <5 4 Ai(£, Ao) + 0(<5 5 ), where Ao £ iR is an explicit eigenvalue 

' 3 For fixed 4 < F 2 < 90, this states that waves are stable for fixed velocity/period and inclination angle 6 sufficiently 
small: equivalently (by rescaling), for fixed inclination angle and period X sufficiently large. 
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associated with the linearized (Bloch) operator for the KdV equation and the corrector Ai(£, Ao) is exactly 
the corrector found in the analogous study of the stability of KdV-KS wavetrains in the singular limit 
5 —> 0 + : see [Barl4, JNRZ15] and Section 1.1.3 above. In particular, there it was proven (through 
numerical evaluation of integrals of certain elliptic functions) that Ind(fc) < 0, as defined in (1.9) for all 
k £ V corresponding to periods X = X(k) in an open interval (X m , Xm) with X m pc 8.44 and Xm ~ 26.1. 
On the other hand, in the regime 0 < | A|/<5 3 +1£| <C 1 a further expansion of the Evans function is needed. 
There, we show that modulo a rescaling of A by <5 3 this expansion is exactly the one derived in [JNRZ15] 
for the singular KdV limit of the KdV-KS equation. From the results of [JNRZ15], this concludes the 
proof of our description of spectral stability in the small-Froude number limit F -A 2 + . Finally, in 
Section 3.2, we carry out a numerical analysis similar to the one in [BJN + 13] where for the KdV-KS 
equation the full set of model parameters was explored: here we consider the influence of 2 < F < oo on 
the range of stability of periodic waves, as parametrized by period X and discharge rate q. 

Authors Note: Since the completion of our analysis, it has been shown that the technical slope condi¬ 
tion (1.4) may indeed be dropped as an assumption to establish the nonlinear modulational stability of 
diffusively stable roll-waves; see [RZ], 


2. Existence and Stability of Roll-Waves in the Limit F —> 2 + 

In this section, we rigorously analyze in the weakly unstable limit F —> 2 + the spectral stability of 
periodic traveling wave solutions of the St. Venant equations (1.10) to small localized (i.e. integrable) 
perturbations. We begin by studying the existence of such solutions and determining their asymptotic 
expansions. In particular, we show that such waves exist and, up to leading order, are described by 
solutions of the KdV-KS equation (1.2) in the singular limit 5 —> 0 (Theorem 1.6). 


2.1. Existence of small-amplitude roll-waves: proof of Thm. 1.6. The goal of this section is to 
establish the result of Theorem 1.6. To begin, notice that traveling wave solutions of the shallow-water 
equations (1.10) with wave speed c are stationary solutions of the system 

(2.1) d t r - d x (u + or) = 0, d t u + d x (j ^2 ~ cu ^j = 1 - ™ 2 + v &x ( r~ 2 d x u) 

of PDE’s. In particular, from the first PDE it follows that u = q — ct for some constant of integration 
q G M and hence r must satisfy the profile ODE 

(2.2) d x [^2 + c2r ) = 1 - r(q - ct) 2 - cvd x ( t~ 2 8 x t) . 


_ j 2 

Clearly then, we have that (■ t,u ) = (to,t 0 7 ) is an equilibrium solution of (2.1) for any to > 0. Further¬ 
more, linearizing the profile ODE (2.2) about r = tq yields, after rearranging, the ODE 


cvtq 2 t" + (c 2 - c 2 ) t’ + 


-- 3/2 


/2-c 


-1/2 


Tn '72 


T = 0, C s = T ° 


- 3/2 


F 


Considering the eigenvalues of the above linearized equation as being indexed by the parameters uq , c, 
and q it is straightforward to check that a Hopf bifurcation occurs when 


c = c q = 


- 3/2 

r o 

F 


and F > 2. 


This verifies that as the Froude number F crosses through F = 2 the equilibrium solutions (r, u) = 
(to,t 0 7 ), corresponding to a parallel flow, becomes linearly unstable through a Hopf bifurcation, and 
hence nontrivial periodic traveling wave solutions of (1.10) exist for F > 2. Moreover, at the bifurcation 
point the limiting period of such waves is given by X = ^ where u = t'^v^^s/F — 2. 

With the above preparation in mind, we want to examine the small-amplitude periodic profiles gener¬ 
ated in the weakly-nonlinear limit F —> 2 + . To this end, we set 5 = \JF — 2 and notice that by rescaling 
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/2, (1-10) 


space and time in the KdV-like fashion Y = 8(x — cq t )/V 1//2 and S = b 3 t/v 1 / 2 , with co = t 0 3 ^ 2 
becomes 

(2.3) 5 2 dsT — dy(u + cor) = 0, b 3 dsu + bdy ~ = ^ 1//2 (l — ™ 2 ) + v x ^ 2 b 2 dy ( r~ 2 dyu ) . 

We now search for small-amplitude solutions of this system of the form (r, u) = (to, t 0 1 ^ 2 ) + S 2 (f, u) with 
wave speed Co in the limit F 2 + . The unknowns f and u satisfy the system 

d 2 dsf — dy(u + cqt ) = 0 

x3o ~ , cq fin o + <5 2 f) 

5 d s u + 6dy { 2pp2 

= v 1/2 b~ 2 (l - (r 0 + b 2 f){r^ X ^ 2 + b 2 u) 2 ^J + u x ^ 2 b 2 dy ((r 0 + b 2 f)~ 2 dyu) . 

Defining the new unknown w = b~ 2 (u + cof) and inserting u = —cqt + b 2 w above yields 



dsT — dyw = 0 


b 3 dsw + 5 x dy 


(r 0 + b 2 f) 2 - Tq + 2t 0 A b 2 f - 3tq b 4 f 2 
2 b 2 F 2 


—3 


+ ( C 0 “ -f2 ) T + 


^5 2 f 2 - 2 c 0 b 2 w 


= V 11 ' 2 ^(2t 0 1/2 C 0 - CoT 0 )f 2 - ‘1t 1 J 2 w - c 0 t 0 2 dyyf + b 2 g(f, w, <5)) 

+ V X / 2 TQ 2 b 2 dyyW + v l l 2 b 2 dy (f(f, b)dy(b 2 VJ — Cof)) 

for some smooth functions g, r. Expanding F~ 2 = | (l — <5 2 ) + 0(<5 4 ) reduces the above system to 


dsf — dyw = 0 

b 3 d s w + bdy ^^3 + ^f 2 - t 0 3/2 w + b 2 f(f, b)^j 

= v 1/2 (^ 2 - 2t 0 1/2 u; - \j^dyyf + b 2 g(f, w, b)\ 

\ 4r ° 2r 0 / 

+ V X / 2 TQ 2 b 2 dy Y W + v l l 2 b 2 dy (f(f, b)dy{b 2 W — Cof)) 


for some smooth function /. Rescaling the independent and dependent variables via 


(Y,S,t,w) (to /4 T, \tq 1/4 5, 3t 0 1 f, 12 tq /2 u>^ , 

we arrive at the rescaled system 


(2.4) 


dsf — dyw = 0 

o - , 

—T7I- ~°SW + 

8t 0 V V/2 


2tq /4 ^ 1 / 2 


dy ( f + \t 2 - w + b 2 f(f , b) 


— \t 2 — W — dyyT + b 2 g(f, w, b) + ^ b 2 dyyW + b 2 dy (r(f, b)dy(b 2 W — CqT )) 


for some smooth functions /, g, and r. 

We now search for periodic traveling waves of the form (f ,w)(Y — cS) in the rescaled system (2.4). 
Changing to the moving coordinate frame (Y — 65,5), in which the 5-derivative becomes zero, and 
integrating the first equation with respect to the new spatial variable Y — 65, we find that w = q — ct 
for some constant q. Substituting this identity into the second equation in (2.4), also expressed in the 
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moving coordinate frame (Y — cS, S), gives 

5 3 5 


yv/2 


c 2 t' + 


8 V v 


2t 0 V V/2 


(1 + c)f + -f 2 + 5 2 /(f,5) 


= —q+ ^t 2 + cf + 5 2 G(t, 5) - ((1 + 5 2 B(t, 5))t)' 

for some smooth functions G and B. Next, introducing the near-identity change of dependent variables 


(2.5) 


T = - f + 6- 


/ o 


B(x, 5)dx 


gives, finally, the reduced, nondimensionalized profile equation 

1, 


( 2 . 6 ) 


T" + -T 2 -cT-q= -5 (c + 1 )T - 


((c+l)T-*r 2 ) + 5 2 m(T,f,5), 5 := 


2t 0 V V/ 2 


for some smooth function m. 

It is well known (see [BD09] for instance) that the limiting 5 = 0 profile equation 

1. 


(2.7) 


To -cT 0 = q 


selects for a given (c, q), up to translation, a one-parameter subfamily of the cnoidal waves of the KdV 
equation ut + uu x + u xxx = 0, which are given given given explicitly as a three-parameter family by 

To(x;ao,k,n) = ao + 12k 2 n 2 cn 2 (ra, k), 

where cn(•,&) is the Jacobi elliptic cosine function with elliptic modulus k G (0,1), k > 0 is a scaling 
parameter, and ao is an arbitrary real constant related to the Galilean invariance of the KdV equation, 
with the parameters (ao, k, k) being constrained by (c, q) through the relations 

c = ao + 4k 2 (2k 2 — 1), 

q = 24fc 2 (l — k 2 )K 4 — ao (^ao + 4k 2 (2k 2 — 1)). 

Note that these cnoidal profiles are 2 I\(k)/n periodic, where K(k) is the complete elliptic integral of the 
first kind. 

Now, noting that (2.6) can be written as 


1 l 


5 \2 


1, 


-(f) 2 + -f 3 --T 2 -qf ) = f ( — ( c + l)f - ^T 2 
6 2 


1, 


+ 6m(f,f,S) 


standard arguments in the study of regular perturbations of planar Hamiltonian systems (see, e.g., [GH90, 
Chapter 4]) imply that, among the above-mentioned one-dimensional family of KdV cnoidal waves To, 
only those satisfying 


( 2 . 8 ) 


r 2K(k)/n 


T' 


c + 1)Tq — — Tq ) dx — 0 


can continue for 0 < 5 <C 1 into a family of periodic solutions of (2.6) and that, further, simple zeros of 
(2.8) do indeed continue for small 5 into a unique, up to translations, three-parameter family of periodic 
solutions of (2.6); that is, for each fixed 0 < 5 <C 1 we find, up to translations, a two-parameter family of 
periodic solutions of (2.6) that may be parametrized by ao and k. The observation that, in the present 
case, the selection principle (2.8) indeed determines a unique wave that is a simple zero, follows directly 
from the proof of Proposition 1.1 (see Remark 1.2) since equation (2.7) implies 

1-2 K{k)/n / 1 \' p2K(k)/n p2 K(k)/n 

J o U Uc + 1)T 0 - -r 2 J dx = J o T' (To + T")' dx = - 1 T ° ( T 0 + To") dx > 

in agreement with the KdV-KS case. This shows that the profile expansion agrees to order 0(5 ) with the 
KdV-KS expansion. Indeed, further computations show that the expansion of the profile coincides with 
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the KdV-KS expansion up to order 0(5 2 ). To complete the proof of Theorem 1.6 we need only observe 
that, instead of fixing the speed as above and letting the period vary, one may alternatively fix the period 
and vary the velocity. 

Remark 2.1. Viewed from a standard dynamical systems point of view, the F —> 2 + limit may be 
recognized as a Bogdanev-Takens , or saddle-node bifurcation; see, e.g., the corresponding bifurcation 
analysis carried out for an artificial viscosity version of Saint Venant in [HC87]. The unfolding of a 
Bogdanev-Takens point proceeds, similarly as above, by rescaling/reduction to a perturbed Hamiltonian 
system [GH90, Section 7.3]. 


2.2. Estimate on possible unstable eigenvalues. Next, we turn to analyzing the spectral stability 
(to localized perturbations) 16 of the asymptotic profiles of the St. Venant equation constructed in Theo¬ 
rem 1.6. To begin, let (fg,ug)(x — cgt) denote a periodic traveling wave solution of the viscous St. Venant 
equation (1.10), as given by Theorem 1.6 for 6 = \/F — 2 E (0, So). More explicitly, in terms of the 
expressions given in Theorem 1.6, the periodic profiles are 


(2.9) 


Tg(0) = To + 5 2l §-TS 


' 5/4 r N 

T » 6 e 


V 


1/2 


/ 5/4,- N 

us(0) = uo + 6 2 fug [\y° 


and the period Xg, wave speed eg, and constant of integration qg = ug + egrg are expressible via 


( 2 . 10 ) 


zz 1 / 2 

^ = ~^ x - 

V 5 


eg — Cq + 


5 2 


4 T, 


3/2 


Cg , 


qg = u 0 + cgT 0 + 


5 2 


12 T ( 


1/2 


q, 


_o /o 

with To, uq constant, Co = r 0 ' /F, and 6 = x — cgt. Linearizing (1.10) about ( fg,ug ) in the co-moving 
coordinate frame 17 ( x — ct , t ) leads to the linear evolution system 

(2.11) dtT — d x (u + cr) = 0, dtu — d x ^cu+ — 2t“ 3 u'^ t^| = — u 2 t — 2 uf + ud x (f~ 2 d x u) 

governing the perturbation (t, u) of (t, u). Seeking time-exponentially dependent modes leads to the 
spectral problem 

(u + ct)' = A T, 

(2.12) 


—3 


v{t 2 u'Y = (A + 2 ut)u - ( ( - 2T 3 Vi' ) t' + cvl ) + ( u 2 - ( - 2t 3 v/ ) ] t. 


T- 3 


F 2 


where primes denote differentiation with respect to x. In particular, notice that (2.12) is an ODE 
spectral problem with V^-periodic coefficients. As described in Section 1.1 above, Floquet theory implies 
that the L 2 (M) spectrum associated with (2.11) is comprised entirely of essential spectrum and can be 
smoothly parametrized by the discrete eigenvalues of the spectral problem (2.12) considered with the 
quasi-periodic boundary conditions ( t,u)(x + Xg) = e^(T,u)(x) for some value of the Bloch parameter 
£ E [—7 t/X g,ir/Xg). The underlying periodic solution (f,u) is said to be (diffusively) spectrally stable 
provided conditions (D1)-(D3) introduced in the introduction hold. Reciprocally, the solution will be 
spectrally unstable if there exists a £ E [— Tr/Xg, ir/Xg) such that the associated Bloch operator has an 
eigenvalue in the open right half plane. 

In this section, we provide a priori estimates on the possible unstable Bloch eigenvalues of the above 
eigenvalue problem (2.12). As a first step, we carefully examine the hyperbolic-parabolic structure of the 
eigenvalue problem and demonstrate that, as F —> 2 + or, equivalently, as 5 —> 0 + , the unstable Bloch 
eigenvalues of this system are 0(1). Next, we prove a simple consistent splitting result that establishes all 
unstable Bloch eigenvalues of (2.12) converge to zero as 5 —> 0 + . We can then bootstrap these estimates to 
perform a more refined analysis of the eigenvalue problem demonstrating that such unstable eigenvalues 
are necessarily 0(5 3 ) as 5 —> 0 + . 


16 The strongest kind of spectral stability, in the sense that it implies spectral stability to co-periodic perturbations, 
subharmonic perturbations, side-band perturbations, etc. 

1 ' Henceforth, we suppress the dependence of f, u q, and c on <5. 
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2.2.1. Unstable eigenvalues converge to zero as 5 —>• 0 + . We begin by ruling out the existence of suffi¬ 
ciently large unstable eigenvalues for (2.12). Setting Z := (t,u,t~ 2 u') t , and recalling that u = q — cf 
for some constant q E M, we first write (2.12) as a first order system 

(2.13) Z'(x) = A(x, A)Z(x), 


where 


(2.14) A(x,X) 


Setting B(x, A) : = 



A/c 

0 

-r 2 /c \ 


0 

0 

/ 2 ) , a 

(q-ci 

■) 2 -cix-a\/c 

A+2 f(q—cf) 

—cr 2 -\-ar 2 /c J 


V 

V 

V / 

A/c 

0 M 



0 

0 r 2 

and noting that A — B is 0(1) as 

a\/c 

V 

V 0 / 




a := t 3 (F 2 + 2 cut'). 


the spectral problem (2.13) is governed by the principal part B(x, A) for |A| sufficiently large. A direct 
inspection shows that the eigenvalues of B are given by 4 and so that the eigenvalues of B 

have two principal growth rates as |A| —> oo. In the following we keep track of both of these spectral 
scales by a series of carefully chosen coordinate transformations preserving periodicity; for details of these 
transformations, see Section 4.1 of [BJRZ11], 

With the above preliminaries, we begin by verifying that the unstable spectra for the system (2.12) are 
0(1) for 5 sufficiently small. Throughout, we use the notation ||u|| 2 = f Q Xs \u(x)\ 2 dx. Note that although 
we focus on uniformity in <5 in the forthcoming estimates of the unstable spectra, the norms || • || do 
depend on <5 through the period X$. 


Lemma 2.2. Let (■ fs,us ) be a family of periodic traveling wave solution of (1.10) defined as in The¬ 
orem 1.6 for all 6 = \JF — 2 £ (0,<5o) for some do > 0 sufficiently small. Then, there exist constants 
Ro,rj > 0 and 0 < <5i < 5q such that, for all 5 £ (0, <5i), the spectral problem (2.12) has no L°°(M) 
eigenvalues with 9i(A) > —p and |A| > Rq. 

Proof. Suppose that A is an L°°(M) eigenvalue for the spectral problem (2.12) and let (t,u) be a cor¬ 
responding eigenfunction satisfying (u, r, t')(Xs) = e^(u, r, r / )(0) for some f £ [—ir/Xg,ir/Xs). As 
described above, setting Z := (t,u,t~ 2 u') t allows us to write (2.12) as the first order system (2.13), 
where the coefficient matrix A(x, A) is given explicitly in (2.14). By performing a series of JW-periodic 
change of variables, carried out in detail in [BJRZ11, Section 4.1], we find there exists a A^-periodic 
change of variables W{-) = P{-\ A, S)Z(-) that transforms the above spectral problem into 

(2.15) W'(x) = ( D x + N)(x, A )W(x), 

supplemented with the boundary condition W(X^) = e 1 ^ W(ff), where the matrices D\,N are defined as 

( 2 . 16 ) D x (; A) = diag^ + 0 o + y, , 9 0 = 

a as in (2.14), and N := ( ^ ) with Np £> a 2 x 2 matrix. Here, 1V(-, A) is an X^-periodic 

V AId,h L!d,d ) 

matrix and, moreover, the individual blocks of the matrix N(-, A) expand as 

Nd,d(’,X) = Np D + A ^-^d^ + A 1 X^ ) D 

(2.17) N D M; A) = K,h + ^K,h 

Nh,d(-,X) = N^ d + A 2 n h,d + A -1 Aj^ D + A 2 Xf 1D , 

with | A r *,* | bounded uniformly in 5 <C 1. Explicit formulae for N J kl and 6\ are given in Appendix A. 
Now, a crucial observation is that, by Theorem 1.6, together with the scalings (2.9)-(2.10), we have 

T 1/2 

lim 9o(x ) = — > 0. 

<5—>-o+ 2 u 
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It follows that do is strictly positive and uniformly bounded away from zero, for all 5 > 0 sufficiently 
small. Thus, there exists 77 , <5i > 0 sufficiently small and Ro > 0 sufficiently large such that if |A| > Ro and 
3ft(A) > — 77 , then the quantity 3 ft(| + $0 + w) is strictly positive and bounded away from zero, uniformly 

in d for 0 < 8 < Si. Likewise, by choosing di smaller if necessary, the quantity |A| 1 / 2 9fi(T'y / /^) may be 
taken to be strictly positive and bounded away from 0 , uniformly in d in the same set of parameters. 
Finally, under the same conditions, taking Ro possibly larger and decomposing 

W:= (W H ,W D ,+,W Dt -) T , 

observing that |W#|, |W / d )+ | and |IFd,-| are .^-periodic functions, it follows by standard energy esti¬ 
mates, taking the real part of the complex L 2 [0, X^j-inner product of each Wj against the ILj-coordinate 
of (2.15)—(2.17), using the above-demonstrated coercivity (nonvanishing real part) of the entries of the 
leading-order diagonal term D\, and rearranging, that there exists a constant C > 0 independent of Ro 
such that, for all 0 < d < Si, 3ft (A) > —77 and |A| > Ro, we have 

\\W H \\ 2 <C(\\W D ,+ \\ + ||!Fd,_||) \\W h \\ and \\W D ,+ \\ 2 + ||1T D ,_|| 2 < CRf 1/2 (||Wb,+ || + ||W A -||) ||W H ||. 

Thus, ||1T//|| < C (||VFd,+ || + ||VFd,-||) < C 2 Rq 1 ' ~\\Wh\\, yielding a contradiction for Ro sufficiently 
large. □ 

Next, we rescale the spatial variable x as y = Sx, noting that, since SXg = z/ 1 // 2 r 0 the period 

is then independent of 5. Then the first-order system (2.13) can be rewritten as 

(2.18) 5Z'(y) = (4>(A) + S 2 A 1 {y-\,8)) Z(y) 

coupled with the boundary condition Z(5Xg) = e^Z(O) for some £ G [— it/SX g, it/SXg), where 

( 2 Atq /2 — 2 Atq /2 0 \ 

^o(A) = j 0 0 Tp zz -1 I 

V A — 2^/to 2A 0 / 

is constant and A, 5) is uniformly bounded (for A and 5 in any compact set). More precisely, for any 
5 in a compact subset of [ 0 , do), we have 

/0(A) 0(A) 0 \ 

ALi(-;A,d)= 0 0 0(1) . 

\ 0 ( 1 ) 0(1) 0 ( 1 )/ 

By analyzing the eigenvalues of ylo(A), we now show that the possible unstable eigenvalues for (2.12) 
converge to the origin as 5 —> 0 + . 

Lemma 2.3. Let ( fg,ug ) be a family of periodic traveling wave solution of (1.10) defined as in Theo¬ 
rem 1.6 for all S = \/F — 2 £ (0, do) f or some do > 0 sufficiently small. Then, for every e > 0, there 
exists a di £ (0, do) such that for all d £ (0, di), the spectral problem (2.12) has no L°°(M) eigenvalues 
with 3?(A) > 0 and |A| > e. 

Proof. By Lemma 2.2, it is sufficient to consider A on a compact set e < |A| < Ro, 3ftA > 0, whence (2.18), 
d —> 0 + represents a uniform family of semiclassical limit problems, with Ao, A\ varying in compact sets. 
By standard WKB-type estimates (see, e.g., the “Tracking Lemma” of [GZ98, ZH98, PZ04]), these have 
no bounded solutions for 0 < d < do sufficiently small, so long as ^4o satisfies consistent splitting, meaning 
that its eigenvalues have nowhere-vanishing real parts: equivalently, j4q(A) has no purely imaginary 
eigenvalue for 3ftA > 0 and e < |A| < Ro- Indeed, it is easy to see using convergence as d —>■ 0 of 
the associated periodic Evans function of Gardner [Gar93] (following from continuous dependence on 
parameters of solutions of ODE), the correspondence between bounded solutions and zeros of the Evans 
function, and analyticity of the Evans function together with properties of limits of analytic functions, 
that A)(A) has a pure imaginary eigenvalue, i.e., the d = 0 version of (2.18) has a bounded solution, if 
and only if there are bounded solutions of (2.18) for a sequence A^ —> A as d —> 0. 
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To prove the lemma therefore, we establish consistent splitting of Ao for A £ A := {A : e < |A|, 9iA > 0}. 
The eigenvalues 7 (A) of the matrix Ao(A) are the solutions of the equation 

(2.19) 7 3 — 2Atq^ 2 7 2 — 2Atq i^ _1 7 + 2A ‘ 1 rJ‘ 2 v~ x + 4Atq = 0. 

Suppose that 7 = iLl £ M?' is an eigenvalue of Aq(A) for some A 6 A. From (2.19) it follows that A must 
be a root of the quadratic equation 

(2.20) 2A 2 Tq^z'~ 1 + 2A ^2tqZz _1 + Q 2 Tq^ — 2iUlTQV~ l ^j — iVL 3 = 0. 

By Lemma 2.2 above, if fl is sufficiently large then the roots of (2.20) satisfy 3t(A) < 0, else, by the 
discussion surrounding the Evans function, above, there would be bounded solutions of (2.18) for SRA > 0, 

| A| large, and h arbitrarily small, a contradiction. (Alternatively, one may repeat the steps of the proof 
of Lemma 2.2 for (2.18) with 5 = 0.) Increasing Q then from the supposed value corresponding to 
an eigenvalue of A\, and tracking the corresponding root A of (2.20), we see that eventually this root 
must cross the imaginary axis in moving from IRA > 0 to SRA < 0. Thus, it is sufficient to search for 
roots of (2.20) of the form A = ?'© for some 0 £ M. Substituting this ansatz into (2.20) and grouping 

o /o 

real and imaginary parts implies that fl and 0 satisfy the system of equations 0 (tq' 0 — Q) = 0 and 
2tq / ' 2 Q = (lx fl 2 /[fl 2 + 2zz _1 Tq //2 ], from which it easily follows that = 0 = 0. It follows that for all 
!1 / 0, the real parts of the roots A j(Q) of (2.20) have constant signs, so that, for each e > 0, Ao(A) 
indeed has consistent splitting in the region A, and the lemma immediately follows. 

□ 

2.2.2. Unstable eigenvalues are 0(h 3 ). Next, we bootstrap the estimates of Lemma 2.2 and Lemma 2.3 to 
provide a second energy estimate on the reduced “slow”, or KdV, block of the spectral problem (2.12) in 
the limit 5 —> 0 + . Notice that this result relies heavily on the fact that the corresponding spectral problem 
for the linearized KdV equation about a cnoidal wave To(-; ao, k, Q{k)) described in Theorem 1.8 has been 
explicitly solved in [BD09, Spe88] using the associated completely integrable structure, and in particular 
has been found to be spectrally stable 18 for all k £ (0,1). 

Proposition 2.4. Let (f$,us) be a family of periodic traveling wave solution of (1.10) defined as in 
Theorem 1.6 for all 5 = \JF — 2 £ (0,ho) f or some 5 o >0 sufficiently small. 

Then there exist positive constants C\, C 2 and S\ £ (0, ho) such that for all 5 £ (0, hi) the spectral problem 
(2.12) has no L°°(M) eigenvalues with R(A) > C\5 A or (5i(A) > 0 and |A| > C 2 S 3 ). 

Proof. By Lemma 2.2 and Lemma 2.3, the possible unstable eigenvalues for the spectral problem (2.12) 
converge to the origin as 5 —X 0 + . To analyze the behavior of these possible unstable eigenvalues further, 
we rescale the unknown Z in (2.18) to W = (Z\, A 2 ^ 3 Z 2 , A 1 / 3 Z 3 ) T , A 1 / 3 denoting the principle third root, 
yielding the system 

SW’(y) = (S 0 (A) + 5 2 A -1 / 3 B\(y \ A, h)) W(y) 

with boundary conditions W(S X$) = e^VF(0) for some £ £ [— tt/SX§, tt/8X§), where 

/2A 2/3 r 0 3/2 — 2 tq /2 0 \ 

Bo(A) = A 1 / 3 I 0 0 Tq zz -1 I 

\A — 2^77 2 A 1 / 3 0 ) 

and A?i(-;A,h) is uniformly bounded for (A ,6) in any compact subset of C x [0, ho). More precisely, for 
any 6 in a compact subset of [0, ho), we have 

/ 0(A 2 / 3 ) 0(A 2 / 3 ) 0 \ 

£i(-;A,h)= 0 0 0(A 2 / 3 ) . 

\0(A 2 / 3 ) 0(1) 0( A 1 / 3 )/ 


18 


’In the Hamiltonian sense, meaning the linearization about To of the KdV equation has purely imaginary spectrum. 
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— 2 T, 


3/2 


M 0 = 


0 

2 


Tn V 


Mi = 


To track the most dangerous terms we write the above system as 

( 2 . 21 ) 5W\y) = (A 1 / 3 M 0 + (2A 2 / 3 + <5 2 X~ 1/3 /3(y)) Mi + \R 0 (X) + 5 2 R^y A, <J)) W(y), 

where 

0 

0 (f 

K -2y/ro 0 0 

and the function /?(•) is some explicit periodic function expressed in terms of To, u and asymptotic 
KdV profiles, while the Rj matrices are uniformly bounded functions of A and 5 on compact subsets of 
C x [0, do)- The goal is to now reduce the first-order problem ( 2 . 21 ) to a constant coefficient problem at 
a sufficiently high order in A and 5. 

To begin, we diagonalize Mq by defining the matrices 



Pn = 


/I 

0 


0 

Kq 


2 7 -; 


3/2 


\ 


0 


0 


uKl 

2t 0 7/ V 



p-1 _ 

M) — 


'1 1 
1 u; 2 

, 1 U! 


0 

V° 


o 


It, 


3/2 


A'o 

0 


0 

0 

2 T, 


\ 


7/2 


_Q_ , 

uK o ) 


where A'o = 4 1 /, 3 Tq/‘V 1 / 3 and u = e 2 * 71- / 3 . Setting Y(-) = P 0 1 1T(-), the system (2.21) can be written in 
equivalent form 

(2.22) 5 Y'{y) = (V/ 3 D 0 - ( 2A 2/3 + 5 2 X~ 1/3 j3(y)) Qi + XR 0 (X) + 8 2 R 1 (y-X,8)\ Y(y ) 

where 

Po = 



\ ( l 

u 

w 2 \ 

j Qi = w 

u? 

1 

1 \cu 2 

1 

0J I 


0 


“A' 0 ur 

and Aj are uniformly bounded in (A, 5) on compact subsets of C x [0,c>o)- This effectively diagonalizes 
system ( 2 . 21 ) to leading order. 

Aiming at reducing (2.22) to a constant-coefficient problem at a higher order, we choose q(-) satisfying 


</(•) = m-m)), 


where (•) denotes average over one period. Notice that the periodicity of /?(•) implies that q is (5 X§)- 
periodic. Now, since the matrix (I 3 + 5 2 X~ l ^q(-)Qi) is invertible for <5 2 A 1//3 sufficiently small, for such 
parameters we can make the change of variables 

U(-) = (l 3 + 5 2 X~ 1 / 3 q(-)Q 1 ) Y(-), 


with U(-) satisfying the first-order system 

5 U'(y) = (A 1 / 3 D 0 - (2A 2 / 3 + <J 2 A’ 1 / 3 (/?(•)» Qi 

+ XR 0 (yX,5 ) +° 5 2 Ri(y, A, 5) + «5 4 A" 2 / 3 R 2 (y A, 5)) U(y), 

with Rj uniformly bounded in (A, 6) on compact subsets of C x [0,5o)- 

Next, we diagonalize the system at a higher order. Since Dq has distinct eigenvalues, we may choose a 
constant matrix Pi such that the commutator [Pi, Dq] equals the off-diagonal part of Q\ . Then provided 
that A and <5 2 A -2 / 3 are small enough we may change the unknown to 

S(-) = (l3+ ^( 2aV3 + 5 2 A" 2 / 3 (/?(•))) Pi) P(-) 

with S'(-) satisfying the first-order system 

(2.23) 8S'{y) = (Pi(A,<5) + XR 0 (yX,8 ) + 5 2 Ri(y, A, 5) + 8 4 A" 1 R 2 (y A, 5)) S(y) 
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where Rj are uniformly bounded in (A, 5) on compact subsets of C x [0, 60 ) and _Di(A, 6 ) is a constant- 
coefficient diagonal matrix whose diagonal entries are 

A*o(A, 5) = —Kq A 1 / 3 - (2A 2 / 3 + 5 2 A’ 1 / 3 (/?(•)}) 




2, ,2 


M+(A , 6 ) = —KqujX 1 / 3 - (2A 2 / 3 + 5 2 A -1 / 3 (/?(•)» 

A- (A, 5) = -I< 0 u 2 A 1 / 3 - ^ (2A 2 / 3 + 5 2 A’ 1 / 3 (/?(■)»• 


In particular, we find 


K(/i 0 (A,<5)) = - (V 0 + ^ <5 2 |A|" 2/3 </?(•)>) ^(A 1/3 ) - 
»(/MM)) = - (V 0 + 3^|- <5 2 |A|- 2 / 3 </?(•)>) K(u;A 1/3 ) - 


2K(AA 1 / 3 ) 


»(/MM)) = - (k 0 + 5 2 |A|“ 2 / 3 </?(•)>) ^A 1 / 3 ) - 

so that, when 9J(A) > 0 and A and 6 2 |A| -2 / 3 are sufficiently small 


2>uKq 

V 

to 

CO 

r 0 

2K(Au;A 1 /3) 

2>vKq 

| A| 2 / 3 ’ 

1 T 0 

2K(Acj 2 A 1 /3) 


ZvKo |A| 2 /3 


®(m+(A ,*))> 
K(M-(A,«5))> 


C|A| 1/3 , 






C'lAI 1 / 3 , 


if 

3(A) 

> 

0 

L ' V |A|2/3 T 

l»(A)|A 
|A| VV ' 

if 

3(A) 

< 

0 ’ 

C'lAI 1 / 3 , 


if 

3(A) 

< 

0 

° V|A|2/3 T 

l»(A)|A 

\M 1/3 J ' 

if 

3(A) 

> 

0 ■ 


Using an energy estimate as in the proof of Lemma 2.2 above, it follows for e, 8 \ >0 sufficiently small 
and R > 0 sufficiently large that for 0 < 5 < system (2.23) has no bounded solutions provided that 

(2.24) »(A)>0, |A| < e, K(A) |A| 1/3 + |3(A)| |A| 2 / 3 > R5 4 . 

In particular, this shows that as <5 —> 0 + the unstable eigenvalues satisfy |A| = 0(5 12//5 ). 

Next, we refine the above bound by using spectral stability of the limiting cnoidal wave. To do so, we 
scale unknowns of the system (2.12) according to 

u{6) = 5 2 ^fb ' 53 


(2.25) 


t( 6) = 5 2 fa 


' 5 / 4 X ^ 


V 


1/2 


' 5/4 ■- > 

5^9 


V 


1/2 


A = 


dTg ^ 4 ^ 1 / 2 


A. 


Our goal is to prove that the rescaled system obtained from (2.12) has no unstable eigenvalues A with 
|A| sufficiently large. To this end, notice that by (2.24) the possible unstable eigenvalues A must satisfy 
the estimate 

(2.26) K(A) > 0, »(A) IAI 1 / 3 + 5 |3(A)| |A| 2 / 3 = 0(1), 

so that, in particular, we already know that A = 0(5 ~ 3 / 5 ). Rewriting the eigenvalue system (2.12) in 
terms of the unknown V(Y) = (5~ 2 (b(Y) + ca(Y)), a(Y), a'(Y)) results in the system 

(2.27) V’(Y) = (A KdV (Y-A) + A5B + 5C(Y) + R(Y;A,5))V(Y) 
with 19 

/ 0 A 0^ 

= I 0 0 1 

\ —1 a o - To(-) 0> 


19 


We don’t need here the exact form of C but we specify it for latter use. 
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B = 


2 t 0 V V/ 2 



C'(-) = 


2tq /4 ^ 1 / 2 


TO 0 

0 0 0 
,0 TO “TO -l + To(-)-a 0y 


and 


/ 0 0 0 

i?(-;A,<S) =10 0 0 

yO(A5 3 ) + 0(5 2 ) 0(A8 3 ) + 0(8 2 ) 0(A8 2 ) + 0(8 2 ) y 

To make use of known results about KdV we rewrite the above spectral problem in a more standard 
way by introducing the unknown IT(-) = (T^(-), V^-), — Vi(-) + (do — To(-)) V 2 {-)) T ■ This leads to 

(2.28) W'(Y) = (H 0 (Y-,A) + A8H 1 +8H 2 (Y) + R(Y-A,8))W(Y) 

with 


(2.29) 


(2.30) H-l = 


and 


0 1 O' 

H 0 (-;A) = ( 0 0 1 

A _ T /(.) ao _ Tq(-) 0, 


1 


^0 0 0 N 
, 1 o 0 
2r 0 i/ V / 2 l 0 o 0, 


0 


#2(0 = 


2tq /4 zv 1 / 2 


TO-TO -i + T 0 (-)-do o 


0 0 0 

i?(-;A,<5) = ( 0(A<5 3 ) + 0(<5 2 ) 0(A8 2 ) + 0(8 2 ) 0(A8 3 ) + 0(8 2 ) 

0 0 0 

The leading order “KdV” part may now be changed to a diagonal constant-coefficient matrix through an 
explicit periodic Floquet change of variable P(-; A) such that 


A 0 0 

P(-;A) = (0 A 1 / 3 0 

,0 0 A 2 / 3 , 


with inverse 


1 


1 1 1 


- 1 u; 2 w + 0(A~ 2 / 3 ) 

3 1 2 


'1 1 1 

1 u co 2 I + 0(A" 2/3 ) 

, 1 CO 2 UJ 

A 0 0 

0 A -1 / 3 0 

,0 0 A -2 / 3 , 


1 U CO" 

Indeed, using that A = 0(8~ 3 / 5 ), replacing IT with Y(-) = P(-; A) _1 IT(-) leads to the system 

(2.31) Y' = (^D(A)+A 2 / 3 8Q 1 +R(-,A,8)'jY 

with 

/ A ^ o n \ 

1 


Vo (A) 0 0 

TO) = | 0 //+( A) 0 

0 0 m-(A), 


1 1 


1 


Q i — 


6 t 0 1 //4 zA/ 2 


, ,2 , ,2 , ,2 

UJ UJ UJ 


UJ UJ UJ 

and R(-',A,8) = 0(8), where the and l-i± are smooth functions of A. In particular, we know that 
when 5R(A) > 0 we have 

»Vo(A)) > 0, »(/i+(A)) < 0, »(/i_(A)) < 0, 


Mo (A) = A 1 / 3 (1 + 0(A- 2 / 3 )) 
M+( A) = A 1 / 3 (a; + 0(A- 2 / 3 )) 
H-(A) = A 1/3 (co 2 + 0(A~ 2/3 )) 
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Next, provided that A 1 / 3 ^ is sufficiently small, which is guaranteed for 8 sufficiently small since A = 
0(6~ 3 / 5 ), we can use a near-identity change of variables to diagonalize the system (2.31) to order 0(A 2 / 3 <5), 
resulting in the system 

Z\Y) = (Di(A,<5) + 0(A 5 / 3 5 4 ) + 0(5)) Z(Y), 
where D i(A,5) is a constant-coefficient diagonal matrix whose diagonal entries are 

Mo(A ,S) = no(A) + ] A 2 / 3 5, 

6t 0 ' v 1 ! 1 

M+(A,5) =/x+(A) + i/4 1/9 A2/3(5 ’ 

6TQ 7 ZV 1 / 2 

MM) = « + (A) + a 2/3 «. 

6t 0 zT/ 2 

Hence, when 3?(A) > 0, A is large, and 5 | A| 1//3 is small, 

^(Mo(A , 6 )) > C'|A| 1 / 3 

»(M+(A,5)) < -O (K(A)|A|- 2 / 3 + 5|A| 2 / 3 ) 

5R(M-(A,5)) < -O (K(A)|A|- 2 / 3 + 8 |A| 2 / 3 ) , 

which, when combined with previous exclusions, is sufficient to prove by an energy argument similar to 
that in Lemma 2.2 above that if <5i G (0, Sq) is sufficiently small and R > 0 is sufficiently large there are 
no eigenvalues of (2.31) when 0 < 8 < 5i and 

M)>o, «tP 5 + |A|2/Sa - R - 

In particular, it follows that any unstable eigenvalue A of the original spectral problem (2.12) must 
satisfy the estimates |A| = 0(5 3 ) and 9i(A) = 0(5 4 ) as 6 -> 0 + , by |A| < i? 3 / 2 and 9i(A) < i?5|A| 2 / 3 < R5, 
together with the scaling A = c5 3 A, with c > 0 a real constant. □ 

2.3. Connection to the KdV-KS index: proof of Thm. 1.8. It follows from Proposition 2.4 that, 
in order to complete the proof of Theorem 1.8, it remains to study the eigenvalues of (2.12), supplemented 
with the appropriate Bloch quasi-periodic boundary conditions, of the form A = A5 3 with A confined 
to a compact subset of C and 0 < 5 <C 1. More precisely, using the rescaling (2.25) from the proof of 
Proposition 2.4 and setting (a,/3) = (5 -2 (6 + ca), a), we must study the spectral problem 

A/3 — a' = 0 

(9 32) 

5 (/3 + t/ 3 + c/3 — a)' = r/3 + c/3 - a - 0" + 0(5 2 )/(a, /3, /?'), 

for some smooth function /, supplemented with the boundary condition (a, /3, f3')(X) = e^(a, f3, (3'){0) 
for some £ G [—7 t/X, tt/X). Here, 8 = \ r 0 1 ^V^ 1 / 2 5 is as in Theorem 1.6, and the bounds in 0(■ ■ ■) are 
uniform as A varies on compact subsets of C. 

To study the above one-parameter family of eigenvalue problems, parametrized by the Bloch fre¬ 
quency £, we can define a periodic Evans function, a complex analytic function whose zeros, for each 
fixed £, agree in location and algebraic multiplicity with the eigenvalues of the boundary value problem 
(2.32). To proceed, we first recall, from the proof of Proposition 2.4, that the spectral problem under 
scaling (2.25) may be written as the dynamical system (Eq. (2.28)) 

W\Y) = ( H 0 {Y-A) + A5H 1 + 5H 2 (Y) + O(5 2 )) W(Y) 

for W = (/ 3,/3',—a + (ctq — To) (3) T , where Hq. Hi, H 2 are given by (2.29) and (2.30). Introducing the 
new dependent variables Z = (W\, W 2 , W 3 + 5((A — Tq — T\)W\ + (—1 + To + cro)H / 2 )) T , the differential 
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system becomes 


(2.33) Z' 


H 0 1 °\ / 0 0 0 \ \ 

0 0 1+5 0 0 0 + 0(5 2 ) Z. 

\ \—A - T' erg — Tq 0/ \T" - T( 2T'-T 1+ A T 0 -l + a 0 ) ) 


Let <!>(•, A, 5) denote the associated 3x3 fundamental solution matrix 20 associated with (2.33). It is an 
easy consequence of the regularity with respect to parameters of the flow associated with the differential 
system (2.33) that <!>(•, A, 5) is analytic with respect to A 6 C and 5. Moreover, for any fixed (£, 5), 
eigenvalues A of (2.32) agree in location and algebraic multiplicity with roots of the Evans function 

(2.34) E sv { A, £, S) := det ($(X, A, 5) - e^ x Id) ; 

see [Gar93] for details. To complete the proof of Theorem 1.8, we must study the roots of the function 
Esv(-,£,5) on compact subsets of C for all £ € [— tt/X, tt/X) and 0 < 5 <C 1. 

In order to connect the stability properties of the small amplitude roll-waves described by Theorem 1.6 
to those of the associated leading order approximating KdV-KS waves given in Proposition 1.1, we also 
define an Evans function ExdV-Ks(E,^, 6 ) for the eigenvalue boundary value problem 

(2.35) A z + ((T s - a 5 )z)'+ z'" + 5 (z" + z"") = 0 

with (z, z 1 , z", z'")(X) = e^ x (z, z', z", z'")( 0). The next result shows that the Evans function ExdV-KS 
is faithfully described, to leading order, by the St. Venant Evans function Egy for 0 < 5 <C 1. 


Proposition 2.5. Uniformly on compact sets of A e C, the Evans function ExdV-KS can be expanded 

for 0 < 5 <C 1 as 

(2.36) 

E Kd v-Ks( A,e,£) = -e* x (l + 0(5))exp(|) (f? sv (A,£, 5) + 0(5 2 (|A| 2 + |£| 2 ) + 0(5 3 (|A| + |£|))) . 

Proof. A similar expansion has been obtained for ExdV-KS i n [JNRZ15, Propositions 3.7 & 4.1]. Our 
proof parallels the arguments there but in a slightly more precise way in order to equate leading order 
terms with those of Egy . 

The starting point is, as in the proof of [JNRZ15, Propositions 3.7], that the spectral problem (2.35) 
may be equivalently written, through a series of variables transformations with Jacobian of size 1 + O(S), 
as 

/ / ° 1 0 \ / 0 0 0 \ \ 

Z' = o 0 1+5 0 0 0 \\z + 0(5 2 ) gi (Z,iu) 

\\-A-T' a 0 -T 0 0/ \T!I - T( 2T' - + A T 0 -l + a 0 J) 


uf = -jjW + 0(d)g 2 (Z, w) 

for some smooth functions gi,g 2 ; see [JNRZ15, Propositions 3.7] for details. To obtain the expected 
homogeneity in (A, £) we also use that the above system supports a conservation law 


AZi 

(2.37) 


[ (To - To - 1 + r 0 + A + 5 (Ti + T[ + T' 1 - T 0 + uq + 1)) Z 1 
+ (To - To - 1 + 5(Ti + 2 T' + A ))Z 2 ) + (-1 - 5(T 0 - t 0 - 1)) Z 3 + 0(5 2 )[Z, w) 


and that the derivative of the profile built in Proposition 1.1 provides a special solution of the system 
when A = 0. Elaborating on this, manipulations on lines and columns of the determinant usually 
performed to validate at the spectral level averaged (or “Whitham”) equations [Ser05] then yield that 
the Evans function EKdV-Ks(E,£, 8 ) may be written as 

(2.38) E KdV - KS ( A, £, 5) = -e l « x (l + 0(5)) exp det [m( A, £, 5) + M(A, £, 5)) , 


20 In particular, this guarantees that 4>(0, A, 5) =Id. 
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where 


M(A,Z,6) = 


N 


\ 0 0 0 


0\ 
0 
0 
1 J 


, N(A,£,8) = 


/ 0(|A| + ICI) 0(1) 0(1) 

0(|A| + |£|) 0(1) 0(1) 

\0(|A| 2 + |A||£|) 0(|A| + |f|) 0(|A| + |e|), 


and 


M( A,£,8) = 


0(5 2 ) 0(5 2 ) 0(5 2 )\ 

0(5 2 ) 0(5 2 ) 0(5 2 

0(<5 2 |A|(|A| + |£|)) 0(5 2 (|A|+ ]£])) 0(h 2 (|A| + ]£])) 0(5 2 

0(5|A|) 0(8) 0(8) 0(8) J 

Again, see [JNRZ15, Propositions 3.7] for details. 

Now, returning to the St. Venant spectral problem (2.33), we notice that (2.33) also supports a con¬ 
servation law in form (2.37) and, furthermore, the space derivative of the traveling-wave profile provides 
a special solution when A = 0. Using calculations completely analogous to those described above for the 
KdV-KS spectral problem (2.35), it follows that we can write 


/ 0(h 2 (|A|+ ]£])) 

o(6 2 m + m 


V 


EsviA,S,8) = (1 + 0(8)) det 

where N is exactly the same matrix as above 21 and 

( 0(<5 2 (|A| + |£|)) 

N(A,U) = 


,U) + N(A,U 


\ 


0(5 2 ) 0(5 2 ) 

0(5 2 (|A| + |e|)) 0(5 2 ) 0(5 2 ) 

VO(5 2 |A|(|A| + |e|)) 0(5 2 (|A| + |e|)) 0(5 2 (\A\ + mj 

The expansion (2.36) follows conveniently by expanding the determinant det(M + M) in (2.38). 


□ 


The proof of Theorem 1.8 now follows immediately from the proof of Proposition 1.3 in [JNRZ15] for 
the stability of KdV-KS waves in the KdV limit 5 —> 0 + . Indeed, the proof in [JNRZ15] followed by 
studying the renormalized KdV-KS Evans function 

E KdV-Ks(A,i,5) = e“^ A (l + 0(h)) exp E KdV - KS (A,£,8 ), 

and using the asymptotic description of E^dV- ks(A,U) up to 0(<5 2 (|A| 2 + lei 2 )) + 0(h 3 (|A| + ]£])). 
Since Proposition 2.5 implies that 

Esv( A, £, S) = E KdV - KS ( A, e, 5) + 0(~5 2 (\ A| 2 + |£| 2 )) + 0(h 3 (| A| + ]£])) 

it follows that the same arguments can be applied without modification to the Evans function Esv( A, £, 5). 
For completeness, we briefly sketch the details. 

For k £ (0,1) such that condition (A) holds, all the non-zero Bloch eigenvalues of the St. Venant 
linearized operator admit a smooth expansion in 5 for 0 < 5 <C 1. In particular, to each pair (£,Ao) 
such that Ao £ Mi is a non-zero eigenvalue of the KdV Bloch operator E^[To] there is a unique root of 
Esv (A, £, 8) for 0 < 8 <C 1 that can be expanded in 5 as 

A(h;e,Ao) = Ao + hA 1 (C,A 0 ) + 0(<5 2 ) 

where Ai(•, •) is the function already involved in definition (1.9) ; see [JNRZ15, Corollary 3.8]. As this 
expansion is uniform in (£, Ao) when (£, Ao) varies in a compact set that does not contain (0, 0), it follows 
that for any neighborhood U of 0 in the spectral plane the condition k £ V provides when 5 is sufficiently 
small a negative upper bound on the real part of the spectrum in C \ {0}; see [JNRZ15, Corollary 3.10] 
for details. 
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Not just in order of magnitude, but by component-wise identification of the coefficients. 


26 



On the other hand, when £ = 0 the origin is an eigenvalue of the KdV Bloch operator £oPo] of algebraic 
multiplicity three and geometric multiplicity two. To unfold the degeneracy for 0 < <5 <C 1 we directly 
apply the arguments in [JNRZ15, Section 4], By studying the three asymptotic regions 0 < 5 < |£| 
[JNRZ15, Lemma 4.4], 5 ~ |£| [JNRZ15, Lemma 4.6], and 0 < |£| < 5 [JNRZ15, Lemma 4.7] separately 
it follows that the condition k G V implies that a set of “subcharacteristic conditions” hold which, in 
turn, are shown to imply that all the roots of Esv{ A,£, 5), with 0 < |(A,£)| <C 1 and 0 < 5 <C 1 all 
have negative real parts. More precisely, under these conditions it is shown that the (diffusive) spectral 
stability conditions (D1)-(D3) and the non-degeneracy hypothesis (HI) are satisfied when 5 is sufficiently 
small. Finally, we note as in the Introduction that, by the description of profiles in Theorem 1.6, slope 
condition (1.4) (or, equivalently, (1.11)) is always satisfied. This completes the proof of Theorem 1.8. 


3. Stability of large amplitude roll-waves 

In the previous section, we rigorously justified the KdV-KS equation (1.2) as a correct description 
for the weak hydrodynamic instability in inclined thin film flow. In particular, in the weakly nonlinear 
regime F —> 2 + we saw that the KdV-KS equation accurately predicts the stability of the associated 
small amplitude roll-wave solutions of the shallow-water equations (1.1). We now complement this study 
by continuing our analysis into the large amplitude regime, far from the weakly-unstable limit F —> 2~*~, 
performing a systematic stability analysis for roll-waves with Froude number on the entire range of 
existence F > 2, including the distinguished limit F —> oo. We begin by considering the limit F -> oo, 
identifying a one-parameter family of limiting systems approachable by various scaling choices in the 
shallow-water equations. We will then numerically study the influence of intermediate Froude numbers 
2 < F < oo on the range of stability of periodic waves. 

3.1. Scaling as F —>■ oo. In our F —> oo studies, we investigate the stability of periodic traveling wave 
solutions of some asymptotic systems obtained from (1.10) via scaling arguments. To begin, performing 
the change of variables (x,t) i-A (x,t) = (k{x — ct),t ) in (1.10) and erasing tildes yields the equivalent 
system 

--2 ' 


dtT — kc d x r — kd x u = 0, dtu — kc d x u + k d x 


2 F 2 


= l — t u z + uk z d x (r 2 d x u). 


To prepare for sending F —> oo above, we first scale the dependent quantities independently via r = F aT a, 
u = F° lu b, k = F° k ko and c = F ac co- This transforms the above shallow-water system to 

F^d t a- k 0 c 0 F° T+ak+ac d x a - k 0 F a - +ak d x b = 0, 

(3.1) F a "d t b- k 0 c 0 F a - +ak+ac d x b + k 0 F Qk ~ 2 ~ 2a - d x 

= 1 - F aT+2au cib 2 + vEl F a -+ 2ak ~ 2a - d x (a~ 2 d x b). 

Seeking stationary 1-periodic solutions of (3.2) leads to an ODE system 
— c 0 F aT+olc a' — F au b' = 0, 

k 0 cl F aT+ak+2ac a + k 0 F ak ~ 2 ~ 2aT = 1 - F a " +2a “a6 2 - i/jfeg c 0 F ac+2ak ~ aT (a^a 1 ) 1 , 

governing the traveling wave solutions of (1.10) in the appropriate moving frame. Note that, by its 
expression as f 0 f(x)dx, the period in the physical Eulerian variables scales as F° T ~ ak under the above 
transformations. To balance the terms of the “elliptic-constant” right-hand side above, we choose now 
to impose a T + 2 a u = 0 and a c + 2ag, — a T = 0. This effectively leaves two free parameters, say a = a T 
and f3 = a with a u = —a/2 and a c = a — 2/3. Under this choice, the above profile system now reads 


- c 0 Fi Q - 2/3 a! - b 1 = 0, 


(3.2) 


feoCg F 3a 3/3 a + ko F 2 2a+/3 1 = 1 — ab 2 — uk 2 } cq (a 2 a)'. 
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The first equation may be integrated to yield b = qo — Fz a ~ 2 P cq a where qo is a constant of integration. 
Substituting this into the second equation above, the rescaled profile system is reduced to the scalar 
equation 

(3.3) k 0 cl F 3a ~ 3 P a! + k 0 F- 2 ~ 2a+ P = 1 - a (q 0 - F^ a ~ 2 ^ c 0 a) 2 - c 0 (a~ 2 a')'. 

In order to ensure that the “elliptic part” ( a~ 2 a')' is not asymptotically negligible as F —> oo, we restrict 
ourselves to the parameter regimes where |a < (3 < 2a + 2, /3 > a. In particular, this latter choice 
requires a > —2. With the above choices, the full rescaled shallow-water equations (3.1) now read as 

F a d t a- koc 0 F 2a ~Pd x a - k 0 F^~ a / 2 d x b = 0, 

(3.4) F- a / 2 d t b-k 0 c 0 F a / 2 -Pd x b + k 0 F~ 2 ~ 2a+ ^ d x 

= 1 — ab 2 + vk$ F 2 P- 5a / 2 d x (a~ 2 d x b ), 

which is a two-parameter family of systems, parametrized by a > —2 and |a < /? < 2a + 2. 

Finally, when sending F —> oo in (3.4), we balance the first-order terms of the reduced profile equation 
(3.2) by requiring 3a — 3/3 = —2 — 2a + /3, i.e. that (3 = \ + 5a/4, effectively reducing our scalings to a 
one parameter family indexed by a > —2. In particular, under this choice the full one-parameter family 
of rescaled shallow-water systems reads as 

F a d t a - k oCo F 3 **/ 4 - 1 / 2 d x a - k 0 F 1/2+3a/4 9 x 6 = 0, 

(3.5) F~ a / 2 dtb — koco i? -3 "/ 4 ” 1 / 2 Q x b + k 0 F” 3 / 2 " 30 / 4 d x 

= 1 — a b 2 + vk$ F d x (a~ 2 d x b), 

and the associated rescaled profile system is equivalent to b = qo — F -1 cq a, qo constant and 

(3.6) F~ 3 / 2 ~ 3a / 4 ^k 0 clci' + k 0 (~Y^J j = 1 - a (q 0 - F~ x c 0 a) 2 - vkl c 0 {a~ 2 a!)'. 

By the discussion above Eq. (3.2), the Eulerian period of the periodic profiles satisfying (3.6) scale as 
p—a/4:— 1/2 ' Notice that (3.6) is equivalent to the profile equation (1.15) claimed in the introduction. 

As mentioned in the introduction, taking F —> oo in (3.6) produces different limiting profile equations 
depending on whether a = —2 or if a > —2. Next, we discuss both limiting profile equations and the 
spectral problems governing the stability of the profiles. 


3.1.1. Case a = —2. Taking a = — 2 in the above discussion corresponds to rescaling (3.2) via 

t = F~ 2 a, u = Fb , k = koF~ 2 , c = cqF 2 

Note that in this case the Eulerian period of the profile is held constant as F —> oo. With this choice, 
the profile equation (3.6) reads as 


(3.7) 


^/cq c 0 a + &0 



1 — a(qo~ F 1 cq a) 2 — vkQ Co (a 2 a')', 


and b = qo~ F 1 Co a, qo constant. Linearizing (3.5) with a = —2 about such an Ao-periodic profile (a, b) 
yields the associated rescaled spectral problem 

F~ 2 \a- koc 0 F~ 2 a'-koF^b' = 0 

FXb — koco F b'—3ko (a -3 a)' 

= —ab 2 — 2 ab + vko F (—2 ah' a + aT 2 b')' 
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(3.8) 



to be considered for (a, b ) satisfying suitable Bloch boundary conditions. Sending F —> oo above, it 
follows that the profile equation is a regular perturbation of the limiting system 


(3.9) 


o! + ko 



i 2 I 2 / —2 /\/ 

1 — a (/o — co (a a ) , 


while the above spectral problem is a regular perturbation of 

Aa — kocoa' — k^o = 0, 

A6 - k 0 c 0 b' - k 0 (a^a)' = -aql + k^v{a^V + 2c 0 a ( 


^oo ®oo 


where (aoo, too); necessarily solutions of (3.9), denote the limiting profiles of (a, 6) as F —»• oo, and b = Fb. 
The limiting profile equation (3.9) is numerically seen to admit periodic orbits existing in a two parameter 
family, parametrized by the period Xq and the discharge rate qo- The stability of these profiles may then 
be investigated by means of the spectral problem (3.10): this is discussed in Section 3.2 below. 


3.1.2. Case a > — 2. When a > —2, it is readily seen that the profile equation (3.6) is a regular pertur¬ 
bation of the ODE 


(3.11) 0 = 1 — q^a — vkQ cq (a 2 a ')', 

which is Hamiltonian in the unknown 1/a. Indeed, denoting h := (// 2 a _1 and rescaling space via x = 
s / \JvkoCoqQ, the above ODE reads as 

(3.12) 0 = 1 — h~ l + h", 

where ’ denotes differentiation with respect to s. This is clearly seen to be Hamiltonian and, upon 
integrating, is equivalent to 

(3.13) n = h - ln(/i) + 

where g is a constant of integration. Elementary phase plane analysis shows that (3.12) admits a one- 
parameter family of periodic orbits parametrized by the constant /./. Indeed, for each n > 1 equation 

(3.13) admits a unique (up to spatial translations) periodic solution whose period we denote X fl . 
Returning to the original variables then, we see that, choosing ko to satisfy vk^coq^X^ = 1, the profile 
equation (3.11) admits a three parameter family of periodic orbits with unit period parametrized by /a, 
co, and qo- Clearly, a necessary condition for such a 1-periodic orbit of (3.11) to persist as a solution of 
(3.6) for F> 1 is that the (a-independent) orthogonality condition 

0 = i (1) X h„cg«' + ko (V)p 1 


be satisfied. Note that this yields a selection principle for the wavespeed via 


Co = 


f (a 1 ) / x (a 2 )'dx [ a 5 x ( a) 2 dx f a 1 x ((a 1 )') 2 dx 

Jo _ _ Jo _ _ Jo _ 

f (a -1 / x adx f a~ 2 x ( a) 2 dx f (( a~ l )') 2 dx 


Jo Jo Jo 

Generically zeros of the above selection function are simple and, in this case, as in Section 2.1 it follows 
from elementary bifurcation analysis that the Hamiltonian profile equation (3.11) admits a two-parameter 
family of periodic orbits, parametrized by /a and qo that persist as periodic orbits of (3.6) for F> 1. In 
particular, note that for a > —2, the specific value of a does not enter into either the limiting profile 
equation nor the selection principle. 

Finally, taking F —> oo in (3.5), we see that the spectral stability of the 1-periodic traveling wave 
solutions (a, qo) of (3.5) constructed above are determined via the spectral problem 


— k 0 b' = 0 , 

= — aqo + vko (2cod~ 3 a a + d~ 2 b')' 
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A a 
0 



where here (a, b ) denotes perturbation, A = F 1 / 2 A, and b = Fb. The former system may also be written 
as 

A a — ko b' = 0 , 

0 = — a cfc + Az//cQCo(a ~ 2 a)' — z/&gCo (a ~ 2 a) w . 

Therefore, in this case we investigate the spectral problem 
(3.14) 0 = h~ 2 a — A a’ + a" . 

for A = (X^)~ l A and d(-) = q q 2 (a -2 a) ((X M ) _1 •) where and X M are associated to a as described 
above. As above, we point out that for a > — 2 , the spectral problem governing the spectral stability 
of the limiting F = oo profiles is independent of the value of a. As a consequence, spectral instability 
of the limiting periodic traveling wave solutions constructed above implies spectral instability of the 
(a-dependent) large-F profiles of the system (3.5) for any a > —2. 

3.2. Numerical investigation as F —> oo. With the above preparations, we report our numerical 
results concerning the existence and spectral stability of the profiles introduced in the previous section 
for the limiting systems in the F —> oo limit. 

In the case a > —2, elementary phase plane analysis indicates that for each /i_ E (0,1) there exists a 
unique (up to translations) periodic solution with h{ 0) = /i_. These profiles were numerically computed 
for 1000 equally spaced values of h- in [0.05,0.95] using the Matlab functions bvp5c and bvp 6 c, with 
absolute and relative error tolerances both set to ICE 8 in the bvp solver. To this end, we utilized a 
bisection method to approximate the value h + > h- such that h+ — \og(h + ) = p and then approximated 

the corresponding period by computing y/2 13 (p — x + log(x ))^ 1 ^ 2 dx. 

As a first attempt to study the L 2 (M)-spectrum associated to (3.14), we utilized a Galerkin truncation 
method known as Hill’s method. For each £ E [— tt/X„, tt/X^), Hill’s method proceeds by expanding 
both the unknown a as well as the background wave h^ in the associated Bloch eigenvalue problem as a 
Fourier series, and then truncating all expansions at some finite order to reduce the problem to finding, 
for each £ E [— tt/X^, ir/X^) the eigenvalues of a finite-dimensional matrix; see Appendix B for more 
details. For each of the profiles numerically constructed above, unstable spectra were present; see Figure 
6 (d) for an example. In addition, we verified the existence of unstable spectra of (3.14) by numerically 
computing winding number for the associated periodic Evans function on a closed contour in the open 
right half complex plane verifying that the winding number is indeed greater than zero for some Bloch 
frequency £. Employing Theorem 1.11, his study suggests that all periodic traveling wave solutions of 
(3.2) are spectrally unstable for F > 1. Note, in our study of spectral problem (3.14) via Hill’s method, 
we used 41 Fourier modes and 1000 Floquet parameters. 



Figure 6 . In (a) and (b) we plot, respectively for F = 38 and F = 100, the spectrum 
corresponding to eigenvalue problem (3.8) with a = —2, and in (c) the spectrum corre¬ 
sponding to the limiting spectral problem (3.10). Here v = 0.1, qo = 0.4, X = 437.6, 
c = 84.0, and (a) Xq = 0.303, (b) Xq = 0.044. In (d), we plot a numerical sampling 
(via Hill’s method) of the unstable spectrum corresponding to the spectral problem 3.14, 
corresponding to the case a > —2, for a representative periodic stationary solution of 
(3.12) . 

Concerning the case a = —2, the profile equation (3.9) was solved by using the Matlab functions bvp5c 
and bvp6c, where we treated co as a free parameter in the bvp solver and used numerical continuation to 
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solve the profile as ko was varied. Our numerical investigation of the associated spectral problem (3.10) 
covered go £ [1.6, 2.2] for v = 2 and go £ [1.2, 2.7] U [0.3,0.45] for v = 0.1 with step size 0.1 in go and 
varying, but smaller, steps in ko in the region where we could solve profiles numerically. Again, we found 
that all these waves have unstable spectra; see Figure 6(c) for an example. In Figure 6, we also plot the 
spectrum as determined using the “a = —2” scaling given in Section 3.1.1 and show that in this scaling 
the spectrum is unstable for large F as well as in the F = oo eigenvalue problem (3.10). Instability in the 
F = oo case was thus confirmed by multiple numerical checks, providing strong evidence of instability. 

In conclusion, our numerical calculations strongly indicate that for F sufficiently large, spectrally stable 
periodic traveling wave solutions of the St. Venant system (1.10) do not exist. This justifies Numerical 
Observation 1 described in 1.2.2. 

3.3. Numerical investigation for intermediate F. In order to better understand the stability of 
periodic traveling wave solutions of (1.10) away from the distinguished limits F —> 2 + and F —> oo, we 
also carried out a numerical investigation for some “intermediate” Froude numbers. 

To solve the profile equation numerically, we expand (2.2) to obtain the profile ODE equation, 

(3.15) t" = (c 2t ' ~ ^3 - 1 + r (« - cr) 2 - 2 czz(t') 2 /t 3 ^ , 

and then proceed the same way as described in the studies of the limiting F —>• oo systems described 
above, this time using a relative error tolerance between 1CD 6 and 1CD 10 in the bvp solver. Checking the 
slope condition (1.11) across numerically determined profiles, we found that it is satisfied for profiles with 
F 3.5 but violated for those with larger F. 

To study the spectrum, we worked with the original (unrescaled) eigenvalue system (2.12), associating 
parameters with those under the rescaling (1.14). We found that the results for Hill’s method were 
more accurate using these original coordinates as opposed to the rescaled coordinates (1.14), even for 
large F. For go = 0.4, ko ranging in the parameter space where profiles could be solved, and for F 
ranging from 10 to either 20 or 30 by 1, we examined, via Hill’s method and the Taylor expansion of the 
Evans function, the cases a = —0.7, —1, —1.4, —1.5, —1.6, —1.7, —1.8, —1.9, —2. We used 201-3000 Fourier 
modes and 21-31 Floquet parameters in Hill’s method and 33-201 Chebyshev nodes for the integral in 
the Taylor expansion. For each value of a, we found that a lower stability boundary curve and an upper 
high-frequency instability curve meet at some value F*(a ) after which no waves are stable. See Figures 4. 
In addition, we find that the upper and lower stability boundaries appear to have a linear relationship of 
the form log {X/u) = b\ log(F) + 62 log(g) + 63; see Figure 7. We used a combination of Hill’s method and 
Evans function computations to determine stability or instability. However, as the period X increases, 
a small loop of spectra parameterized by the whole interval of Floquet parameters £ £ [— tt/X, n/X) 
shrinks until eventually neither the Evans function nor Hill’s method can definitively determine small 
frequency stability at which point we rely on Hill’s method for the overall behavior of the spectrum; in 
this region, analytic verification of the stability regions as F —» 00 would be beneficial. However, for the 
lower intermediate F region where the period is smaller, the stability picture appears clear. In particular, 
for a = —2 our numerical study, though not a numerical proof, strongly suggest the stability region 
shown in Figure 4 (a). As F continues to increase, it becomes easier again to compute the spectrum 
with Hill’s method as the spectrum approaches that of the limiting system given in equation (3.9); see 
Figure 6 (a)-(c) for an illustration. For F = 38, we examined the full set of periods corresponding to 
those examined in Figure 4, and found no stability region, confirming that by this point the upper and 
lower stability boundaries have met. 

Acknowledgement: The numerical stability computations in this paper were carried out using the 
STABLAB package developed by Jeffrey Humpherys and the first and last authors; see [BHZ] for docu¬ 
mentation. We gratefully acknowledge his contribution. We thank also Indiana University Information 
Technology Service for the use of the Quarry computer with which some of our computations were car¬ 
ried out. Authors in various combinations also acknowledge the kind hospitality of ENS Paris, Indiana 
University, INSA Toulouse and Universite Lyon 1 that hosted them during parts of the preparation of 
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Figure 7. In (a) and (b), black dots mark the computed boundary and pale dots (red 
in color plates) mark the best least curve fit. (a) Lower stability boundary. We have 
log(X) «= b\ log(F) + £> 2 log(< 7 ) + 63 where b\ = —0.692, 62 = 3.46, and 63 = 0.3. Here 
a = —1.6, —1.8, —1.9, —2. The maximum error is 0.200 and the maximum relative error 
is 0.056. The average relative error is 0.012 and the average absolute error is 0.041. (b) 
High frequency stability boundary. We have log(X) ~ 61 log(F) + & 2 log(g) + 63 where 
b\ = —0.791, 62 = 1.73, and 63 = 3.92. Here a = —1.6,—1.8,—1.9,—2. The maximum 
error is 0.228 and the maximum relative error is 0.052. The average relative error is 0.024 
and the average absolute error is 0.103. 
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Figure 8. Plot of spectra for the intermediate eigenvalue problem. We plot the ap¬ 
proximations returned by Hill’s method with black dots and those obtained by Taylor 
expanding the Evans function with light dashes. We plot the spectra corresponding to 
Floquet parameter zero with a pale star (red in color plates). The solid line indicates the 
imaginary axis. In all cases, a = —2, v = 0.1, and qo = 0.4. Other parameters are (a) 
F = 6, X = 7.83, (b) F = 6, X = 8.78, (c) F = 10, X = 76.9, (d) F = 10, X = 90.9. 
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Appendix A. Notation for High Frequency Bounds 


The proof of Lemma 2.2 in Section 2.2.1 relied heavily on computations previously carried out in 
detail in Section 4.1 of [BJRZ11]. There, the authors were concerned with the stability analysis of 
traveling solitary wave solutions of the viscous St. Venant equation (1.10), i.e. those traveling wave 
solutions that decay exponentially fast to zero as x — cf -> Too. By a straightforward adaptation 
of this analysis to the periodic case, 22 it follows that, for each k € V, there exists an A^-periodic 
change of variables W(-) = P(-; A, S)Z(-) that transforms the spectral problem (2.13) into an equiv¬ 
alent system of the form W'{x) = (D(x, A) + 0(x, A)) W(x), supplemented with the boundary condi¬ 
tions W(Xs) = e^lT(O) for some £ G [— it/Xs, n/Xs), where the 3x3 matrix-valued D is defined 

via £>(•, A) = diag + 9 0 + 


™* 2 i ™ _ r 2 (q— ct) 2 


where 9q = X— and 9i = - 1 


a := t S (F 2 + 2 cut'). Moreover, the 3x3 matrix 0(x, A) has the block structure 0 = 


©+- 

0 _+ 0 _ 

where 0++ is a 2 x 2 matrix. The individual blocks of the matrix 0 can be expanded as cubic polynomials 
in A -1 / 2 with matrix valued coefficients. More precisely, they expand as 


— -^^72 with 

yJV V 6 ' Z 


©++ 


©++(-, A) = 0° + + A _1 / 2 0^ + + A' 1 © 2 + + A" 3 / 2 © 3 + 
©+_(-, A) = 0° _ + A“ 1/2 0^_ + A" 1 © 2 _ + A" 3 / 2 © 3 _ 
©_+(-, A) = 0° + + A“ 1/2 ©L + + A" 1 © 2 + + A" 3/2 0 3 + 
[ ©__(-, A) = 0° _ + \- 1/2 Q l __ + A” 1 © 2 _ + A“ 3/2 0 3 _, 


with 


0 ° - 
u ++ — 


©++ = 


0 2 - 


0 3 — 
u ++ — 


©+- = 


0i_ = 


2v r 2 v J 
0 

i/Z ((g - ct ) 2 + ca^y 

o 2 f 3 {q-cf) 

_ v 

f 2 (q—cf) 


,0 


2 v 1 

f 3 (q—cf) 2 

T 3 / 2 

0 

2f 3 (q—cf) 


COL r 

2^/v v 3 ! 2 , 


v 2 


T 2 (q-cr ) 2 . ca f 3 

^2^v 3 ! 2 , 


2v 


_ / OLT 2 T CT 2 

_ ' ~2v ~ 


ct• 
2v 


0 : 


0 L = (A 


-\ 2 |- T 2 

q - ct) + ca — 


Q\ = (_ r 2 (g-cf) + ca 


_ 

2y/u z/ 3 / 2 


^_2t_ I r 3 /a i -\ 

VL + v 3 / 2 V c + c > 
f 2 {q—cf) _ a f 3 

\ 2 v 3 / 2 


"^ + X§+* 

T 2 (q—CT ) . ce f 3 
JV ' 2 z/3/2 


©i- = 


©i_ = 


f 2 

V 


\2f 2 v 

f 3 (q—cf) 2 


1 / 3/2 


0 


sjv 1 


-=(° 
r 2 (q—cr ) 


T 2 (g—cr) 2 COL 


/v 

*2 


2v 

a f 3 
2 z/ 3/2 


f 3 

2/? I/ 3 / 2 

0 ° _ = 


C)2: = ( T 2 (q-cr) 2 

- ' 2i> 


ca t 3 
2 yd/ 1 / 3/2 


2z/ 


From this, the matrices Ny> t D, Nh,d-, Xd,h hi (2.17) are obtained through the identification 


©++ 

©-+ 


©H— 

0 _ 


0 Nh,d 
Xd,h X D j) 


where A^dd is a 2 x 2 matrix. 


22 


Namely, using formally identical coordinate changes depending on the profile and its derivatives. 
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Appendix B. Computational Methods 


For completeness, we very briefly describe the computational methods utilized in our investigations 
reported in Section 3.2 and 3.3 above. For more details, the interested reader is referred to [BJN+13] 
where an analogous numerical study is performed on the generalized Kuramoto-Sivashinsky equation. 

B.l. Hill’s method. To determine the global picture of spectrum of a linear A-periodic operator L, 
we use Hill’s method. The linear operator L takes the form Lj^ fj,k,q(x where the fj : k,q(x ) 

are A periodic. Following [DKCK07], we represent the coefficient functions fj t k,q( x ) as a Fourier series 
fj,k,q(x ) = YIJL-oo ^j,k,q el2n ^ x ^ X ■ We use Matlab’s fast Fourier Transform to determine the coefficients 
4>j,k,q- The generalized eigenfunctions are represented as 23 v(x ) = e *? x v ] e mjx l x , where £ £ 

[— ir/2X, n/2X) is the Floquet exponent. Upon substituting the Fourier expansions into the eigenvalue 
problem, fixing £, and equating the coefficients of the resulting Fourier series, we arrive at the eigenvalue 
problem L^v = Xv where is an infinite dimensional matrix. The spectrum of the operator L is given 
by cr(L) = U ?e[ _ -n/ 2 X tt/ 2 X) Truncating the Fourier series at N terms leads to a finite dimensional 

eigenvalue problem L^ N v = Xv. The matrix L^ N is of the form 1 M i where M\v = XM 2 V is the 
original eigenvalue problem. Typically M 2 is the identity, but in (3.14), M 2 is diagonal with jth diagonal 
entry i(j + £), hence we avoid £ = 0 in that case so that M 2 is invertible. We compute (t(M n ) on a 
mesh to approximate the spectral curves of L. For our numerical studies, we used the implementation 
of Hill’s method built into STABLAB [BHZ], For discussion of Hill’s method and its convergence, see 
[CD10, DK06, JZ12]. 

B.2. Evans function. Our results for Hill’s method are augmented by use of the Evans function. To 
this end, note all the spectral problems we study, such as theone given in (3.14), may be written as a 
first order system of the form W’{x) = A(x; A)IU(x) and that, further, A £ C belongs to the essential 
spectrum of the associated linearized operator L if and only if this 1st order system admits a non-trivial 
solution satisfying 

Y(x + A; A) = e^ x Y{x; A), Vx £ R 

where here A denotes the period of the coefficients of L. Following Gardner [Gar93], the Evans function 
is defined as D{ A,£) := det (T(A, A) — e l ^ x ) where the matrix ^(a^A) satisfies d x ^(x, A) = Af(i, A) 
and ^(CfA) = Id. By construction then, the roots of D (•;£) agree in location and algebraic multiplicity 
with the eigenvalues of the associated ^-dependent spectral problem. Unfortunately, the Evans function 
as described here is poorly conditioned for numerical computation. To remedy this, as in [BJN + 13], we 
use the observation of Gardner [Gar93] that 

D( A, 0 := det(T(A) - e^'ld) = det e '^ Id ) , 

to express the Evans function as an exterior product of solutions of 

with data (Id,Id) T at x = 0 and (e^Id,Id) T at x = A; for details see [BJN+13]. We then use the 
polar coordinate method of [HZ06] to evolve the solutions. This algorithm is numerically well-conditioned 
[Zum09]. All computations were carried out using STABLAB [BHZ], 

As mentioned above, Hill’s method is ideal for obtaining a global picture of the spectrum and the Evans 
function can be evaluated on contours and the winding number evaluated to determine the presence of 
zeros. However, neither method determines definitively whether unstable spectra of arbitrarily small size 
exist, due to numerical error. In particular, such methods can not be used to determine the spectrum of 
the associated linearized operators in a sufficiently small neighborhood of the origin, i.e. they can not 

2 ^Mark that in the standard implementation of Hill’s method a periodic wave is treated as a periodic function of twice 
its fundamental period. As recalled in [Rodl3, Section 3.1, p.67], this is originally motivated by the fact that in applications 
to self-adjoint second-order scalar operators, the Floquet-zero spectrum will then provide edges of spectral bands. 
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resolve the modulational instability problem. A strategy for rigorously computing stability , which we 
utilize in our intermediate F numerical studies reported in 3.3 above, involves a Taylor expansion of the 
Evans function as we now briefly describe; see [BJN + 13] for details. 

Due to the presence of a conservation law in the governing system (see [Ser05, JZN11, JNRZ14, Rodl3]) 
the Evans function has a double root at the origin when £ = 0. As such, the Taylor expansion of the 
Evans function about the origin (A, £) = (0, 0) takes the form D( A, £) = C 2 ,oA 2 + cipA^ + co^C 2 + C 3 ,oA 3 + 
C 2 ,iA 2 £ + ci^A^ 2 + co, 3 ^ 3 + 0(|A| 4 + |£| 4 ) where the coefficients Ckj may be determined via Cauchy’s 
integral formula, 

(B.l) Cfcj . = -_L/ / i>(A,0A- fc - 1 r J '“ 1 rfAde 

47r JdB(0,r) J8B(0,r) 


with r > 0 sufficiently small. Setting otj = 


—Cl,l + (— iy +1 yjc\ t - 4C2,0C0,2 


> Pj ~ 


C3,0O?+C2,ia^+Cl,2aj+C0,3 


3 2c 2 ,o ’ ^3 2c2,ootj+ci,i 

readily checks that the roots of the Evans function near (A, £) = (0,0) may be continued for |£| < 1 as 


one 


A j(0 — aji + /3j^ 2 + f (1 — s)"\”'(s^)ds. 
1 Jo 


The spectral curves at the origin are thus approximated by ay£+/3j£ 2 with spectral stability corresponding 
to the case otj G Mi and < 0; see [BJN + 13] for details. 

In practice, to compute the Taylor expansion coefficients, rather than compute the Evans function on 
the contour integral in the variable A for fixed £ given in (B.l), we interpolate with J2k—o F k ^ x (K = 3 is 
the largest power of e^ x that appears) and then use the Taylor expansion e lk ^ x = l + (iA’£.T) + (iA;£x) 2 /2+... 
yielding D( A,£) = YlkLo c k£ k i from which the contour integral can be determined simply by reading off 
the corresponding coefficient. Calling the quantity just determined D, we see 


1 

27u 


Dj A) 

\\\=R 1 A r+1 


dX = 


1 

27T 


i9 \ 


D(Re 

R r e ire 


dd = 


2 R r 


D(Re™ u )e~ ir7rt) d6, 


i -1 


which we compute by approximating the integrand with Chebyshev interpolation and integrating. 


Appendix C. Computational effort 

C.l. Computational environment. In carrying out our numerical investigations, we used a Mac.Book 
laptop with 2GB memory and a duo core Intel processor with 2GHz processing speed, a 2009 Mac Pro 
with 16GB memory and two quad-core intel processors with 2.26 GHz processing speed, and Quarry, a 
supercomputer at Indiana University consisting of 140 IBM HS21 Blade servers with two 2GH quad-core 
Intel Zeon 5335 processors per node and delivering 8.96 teraflops processing speed. All computations 
were done using Matlab and the Matlab based stability package STABLAB. 

C.2. Computational time. We begin by providing computational statistics for the representative pa¬ 
rameter set a = —2, qo = 0.4, F = 10, and X = 50. We compute the Evans function, D( A, £), on a semi¬ 
circular contour, H, of radius R = 0.2 with 42 evenly spaced Floquet parameters £ £ [—7 t/X, — 7r/(10AT)] U 
[7 t/(10AT), toIX]. We require the relative error between contour points of the image contour, Re, D(-,£) : 
II —> Y(:, not exceed 0.2 so that Rouche’s theorem implies the winding number of corresponds to the 
number of roots of D(-,£) in H. We use 277 points in D, chosen adaptively, to achieve 0.2 relative error 
in each Yg at a computational cost of 56.0 seconds using 8 Matlab workers on Quarry to determine the 
winding number is zero. Computing the Taylor expansion of the Evans function at the origin requires 
61.7 seconds on Quarry, while computing the spectrum via Hill’s method using 603 Fourier modes and 
21 Floquet parameters comes at a computational cost of 143 seconds. 

As the period X increases or as F increases, the number of Fourier modes needed in Hill’s method 
increases. Using 600 Fourier modes typically takes around 3 minutes on the Mac Pro, while using 1600 
Fourier modes takes about 77 minutes, and using 3000 Fourier modes requires approximately 8 hours. 

In creating Figure 4 (a) it took 4.36 days of computation time to evaluate the Taylor coefficients and 
34.5 days to compute the spectrum using Hill’s method, while the Evans function required an estimated 
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58 hours. A typical profile requires only a few seconds to compute, but we must use continuation whereby 
we use a nearby profile as an initial guess in the boundary value solver, so that computing the profiles 
also required a great computational effort. Overall, taking into account the use of parallel computing and 
all values of a investigated, we estimate that total computational time for the project exceeds a year. 


[AeM91] 

[Barl4] 

[BD09] 

[BGMR] 

[BGNR14] 

[BHZ] 

[BJN+10] 

[BJN+12] 

[BJN+13] 

[BJRZ11] 

[BL05] 

[BM04] 

[BN95] 

[Bro69] 

[Bro70] 

[CD02] 

[CD10] 

[CDK93] 

[CE90] 

[DK06] 

[DKCK07] 

[Dre49] 

[EMR93] 

[FSMA03] 

[FST86] 

[Gar93] 

[Gar97] 


References 

Mina B. Abd-el Malek. Approximate solution of gravity-affected flow from planar sluice gate at high froude 
number. J. Comp. & Appl. Math., 35(13):83 - 97, 1991. 

Blake Barker. Numerical proof of stability of roll waves in the small-amplitude limit for inclined thin film flow. 
J. Differential Equations, 257(8):2950-2983, 2014. 

Nate Bottman and Bernard Deconinck. KdV cnoidal waves are spectrally stable. Discrete Contin. Dyn. Syst., 
25(4): 1163—1180, 2009. 

Sylvie Benzoni-Gavage, Colin Mietka, and L. Miguel Rodrigues. Co-periodic stability of periodic waves in some 
Hamiltonian PDEs. In preparation. 

Sylvie Benzoni-Gavage, Pascal Noble, and L. Miguel Rodrigues. Slow modulations of periodic waves in Hamil¬ 
tonian PDEs, with application to capillary fluids. J. Nonlinear Sci., 24(4):711-768, 2014. 

Blake Barker, Jeffrey Humpherys, and Kevin Zumbrun. STABLAB: A MATLAB-based numerical library for 
Evans function computation. Available at: http://impact.byu.edu/stablab/. 

Blake Barker, Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, and Kevin Zumbrun. Whitham averaged 
equations and modulational stability of periodic traveling waves of a hyperbolic-parabolic balance law. Journes 
Equations aux derivees partielles, pages 1-24, 6 2010. Available as http://eudml.org/doc/116384. 

Blake Barker, Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, and Kevin Zumbrun. Stability of periodic 
Kuramoto-Sivashinsky waves. Appl. Math. Lett., 25(5):824-829, 2012. 

Blake Barker, Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, and Kevin Zumbrun. Nonlinear modula¬ 
tional stability of periodic traveling-wave solutions of the generalized Kuramoto-Sivashinsky equation. Phys. D, 
258(0):11 - 46, 2013. 

Blake Barker, Mathew A. Johnson, L. Miguel Rodrigues, and Kevin Zumbrun. Metastability of solitary roll wave 
solutions of the St. Venant equations with viscosity. Phys. D, 240(16):1289-1310, 2011. 

Abdellatif Boudlal and V. Yu Liapidevskii. Stability of regular roll waves. J. Computational Technologies, 10 
(2):3-14, 2005. 

Neil J. Balmforth and Shreyas D. Mandre. Dynamics of roll waves. J. Fluid Mech., 514:1-33, 2004. 

Doron E. Bar and Alexander A. Neponmyashchy. Stability of periodic waves governed by the modified Kawahara 
equation. Phys. D, 86(4):586-602, 1995. 

Robert R. Brock. Development of roll-wave trains in open channels. J. Hydraul. Div., 95(4):1401-1428, 1969. 
Robert R. Brock. Periodic permanent roll waves. J. Hydraul. Div., 96(12):2565-2580, 1970. 

Hsueh-Chia Chang and Evgeny A. Demekhin. Complex wave dynamics on thin films, volume 14 of Studies in 
Interface Science. Elsevier Science B.V., Amsterdam, 2002. 

Christopher W. Curtis and Bernard Deconinck. On the convergence of Hill’s method. Math. Comp., 79(269):169- 
187, 2010. 

Hsueh-Chia Chang, Evgeny A. Demekhin, and Dmitry I. Kopelevich. Laminarizing effects of dispersion in an 
active-dissipative nonlinear medium. Physical Review D, 63(3-4):299-320, 1993. 

Pierre Collet and Jean-Pierre Eckmann. Instabilities and fronts in extended systems. Princeton Series in Physics. 
Princeton University Press, Princeton, NJ, 1990. 

Bernard Deconinck and J. Nathan Kutz. Computing spectra of linear operators using the Floquet-Fourier-Hill 
method. J. Comput. Phys., 219(1):296-321, 2006. 

Bernard Deconinck, Firat Kiyak, John D. Carter, and J. Nathan Kutz. SpectrUW: a laboratory for the numerical 
exploration of spectra of linear operators. Math. Comput. Simulation, 74(4-5):370-378, 2007. 

Robert F. Dressier. Mathematical solution of the problem of roll-waves in inclined open channels. Comm. Pure 
Appl. Math., 2:149-194, 1949. 

Nicholas M. Ercolani, David W. McLaughlin, and Heinz Roitner. Attractors and transients for a perturbed 
periodic KdV equation: a nonlinear spectral analysis. J. Nonlinear Sci., 3(4):477-539, 1993. 

B. Freeze, Sergey Smolentsev, Neil Morley, and Mohamed A. Abdou. Characterization of the effect of froude 
number on surface waves and heat transfer in inclined turbulent open channel water flows. Int. J. Heat & Mass 
Transfer, 46(20):3765 - 3775, 2003. 

Uriel Frisch, Zhen-Su She, and Olivier Thual. Viscoelastic behaviour of cellular solutions to the Kuramoto- 
Sivashinsky model. J. Fluid Mech., 168:221-240, 1986. 

Robert A. Gardner. On the structure of the spectra of periodic travelling waves. J. Math. Pures Appl. (9), 
72(5):415-439, 1993. 

Robert A. Gardner. Spectral analysis of long wavelength periodic waves and applications. J. Reine Angew. Math., 
491:149-181, 1997. 


36 



[GH90] 

[GZ98] 

[HC87] 

[HZ06] 

[Jef25] 

[JNRZ14] 

[JNRZ15] 

[JZ10] 

[JZ11] 

[JZ12] 

[JZN11] 

[KSF84] 

[KT75] 

[Kur84] 

[Mie97a] 

[Mie97b] 

[Mie02] 

[Nob06] 

[NR13] 

[OZIO] 

[PSU07] 

[PZ04] 

[RG12] 

[RG13] 

[Rodl3] 

[RZ] 

[Sch96] 

[Sch98] 

[Ser05] 


John Guckenheimer and Philip Holmes. Nonlinear oscillations, dynamical systems, and bifurcations of vector 
fields, volume 42 of Applied Mathematical Sciences. Springer-Verlag, New York, 1990. Revised and corrected 
reprint of the 1983 original. 

Robert A. Gardner and Kevin Zumbrun. The gap lemma and geometric criteria for instability of viscous shock 
profiles. Comm. Pure Appl. Math., 51(7):797-855, 1998. 

Shyh Hong Hwang and Hsueh-Chia Chang. Turbulent and inertial roll waves in inclined Him flow. Phys. Fluids, 
30(5): 1259-1268, 1987. 

Jeffrey Humpherys and Kevin Zumbrun. An efficient shooting algorithm for Evans function calculations in large 
systems. Phys. D, 220(2):116-126, 2006. 

Harold Jeffreys. Lxxxiv. the flow of water in an inclined channel of rectangular section. Philos. Mag. Ser. 6, 
49(293):793-807, 1925. 

Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, and Kevin Zumbrun. Behavior of periodic solutions of 
viscous conservation laws under localized and nonlocalized perturbations. Invent. Math., 197(1): 115—213, 2014. 
Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, and Kevin Zumbrun. Spectral stability of periodic wave 
trains of the Korteweg-de Vries/Kuramoto-Sivashinsky equation in the Korteweg-de Vries limit. Trans. Amer. 
Math. Soc., 367(3):2159-2212, 2015. 

Mathew A. Johnson and Kevin Zumbrun. Nonlinear stability of periodic traveling wave solutions of systems of 
viscous conservation laws in the generic case. J. Differential Equations, 249(5):1213-1240, 2010. 

Mathew A. Johnson and Kevin Zumbrun. Nonlinear stability of periodic traveling-wave solutions of viscous 
conservation laws in dimensions one and two. SIAM J. Appl. Dyn. Syst., 10(1): 189—211, 2011. 

Mathew A. Johnson and Kevin Zumbrun. Convergence of Hill’s method for nonselfadjoint operators. SIAM J. 
Numer. Anal, 50(l):64-78, 2012. 

Mathew A. Johnson, Kevin Zumbrun, and Pascal Noble. Nonlinear stability of viscous roll waves. SIAM J. Math. 
Anal., 43(2):577-611, 2011. 

Evegenii A. Kuznetsov, M. D. Spector, and Gregory E. Fal'kovich. On the stability of nonlinear waves in integrable 
models. Phys. D, 10(3):379-386, 1984. 

Yoshiki Kuramoto and Toshio Tsuzuki. On the formation of dissipative structures in reaction-diffusion systems, 
reductive perturbation approach. Prog. Theor. Phys., 54(3):687-699, 1975. 

Yoshiki Kuramoto. Chemical oscillations, waves, and turbulence, volume 19 of Springer Series in Synergetics. 
Springer-Verlag, Berlin, 1984. 

Alexander Mielke. Instability and stability of rolls in the Swift-Hohenberg equation. Comm. Math. Phys., 
189(3):829-853, 1997. 

Alexander Mielke. Mathematical analysis of sideband instabilities with application to Rayleigh-Benard convec¬ 
tion. J. Nonlinear Sci., 7(l):57-99, 1997. 

Alexander Mielke. The Ginzburg-Landau equation in its role as a modulation equation. In Handbook of dynamical 
systems, Vol. 2, pages 759-834. North-Holland, Amsterdam, 2002. 

Pascal Noble. On the spectral stability of roll-waves. Indiana Univ. Math. J., 55(2):795-848, 2006. 

Pascal Noble and L. Miguel Rodrigues. Whitham’s modulation equations and stability of periodic wave solutions 
of the Korteweg-de Vries-Kuramoto-Sivashinsky equation. Indiana Univ. Math. J., 62(3):753-783, 2013. 
Myunghyun Oh and Kevin Zumbrun. Stability and asymptotic behavior of periodic traveling wave solutions of 
viscous conservation laws in several dimensions. Arch. Ration. Mech. Anal., 196(1): 1 —20, 2010. Erratum: Arch. 
Ration. Mech. Anal., 196(l):21-23, 2010. 

Robert L. Pego, Guido Schneider, and Hannes Uecker. Long-time persistence of Korteweg-de Vries solitons as 
transient dynamics in a model of inclined Him flow. Proc. Roy. Soc. Edinburgh Sect. A, 137(1):133-146, 2007. 
Ramon Plaza and Kevin Zumbrun. An Evans function approach to spectral stability of small-amplitude shock 
profiles. Discrete Contin. Dyn. Syst., 10(4):885-924, 2004. 

Gael L. Richard and Sergey L. Gavrilyuk. A new model of roll waves: comparison with Brock’s experiments. J. 
Fluid Mech., 698:374-405, 2012. 

Gael L. Richard and Sergey L. Gavrilyuk. The classical hydraulic jump in a model of shear shallow-water flows. 
J. Fluid Mech., 725:492-521, 2013. 

L. Miguel Rodrigues. Asymptotic stability and modulation of periodic wavetrains, general theory & applications 
to thin film flows. Habilitation a diriger des recherches, Universite Lyon 1, 2013. 

L. Miguel Rodrigues and Kevin Zumbrun. Periodic-coefficient damping estimates, and stability of large-amplitude 
roll waves in inclined thin film flow. In preparation. 

Guido Schneider. Diffusive stability of spatial periodic solutions of the Swift-Hohenberg equation. Comm. Math. 
Phys., 178(3):679-702, 1996. 

Guido Schneider. Nonlinear diffusive stability of spatially periodic solutions—abstract theorem and higher space 
dimensions. In Proceedings of the International Conference on Asymptotics in Nonlinear Diffusive Systems 
(Sendai, 1997), volume 8 of Tohoku Math. Publ., pages 159-167, Sendai, 1998. Tohoku Univ. 

Denis Serre. Spectral stability of periodic solutions of viscous conservation laws: large wavelength analysis. 
Comm. Partial Differential Equations, 30(l-3):259-282, 2005. 


37 



[Siv77] 

[Siv83] 

[Spe88] 

[SS01] 

[Win93] 

[YY 03] 

[ZH98] 

[Zum07] 

[Zum09] 


Gregory I. Sivashinsky. Nonlinear analysis of hydrodynamic instability in laminar flames. I. Derivation of basic 
equations. Acta Astronaut., 4(11-12):1177-1206, 1977. 

Gregory I. Sivashinsky. Instabilities, pattern formation, and turbulence in flames. Ann. Rev. Fluid Mech., 
15(1):179-199, 1983. 

M. D. Spektor. Stability of conoidal [cnoidal] waves in media with positive and negative dispersion. Zh. Eksper. 
Teoret. Fiz., 94(l):186-202, 1988. 

Bjorn Sandstede and Arnd Scheel. On the stability of periodic travelling waves with large spatial period. J. 
Differential Equations, 172(1): 134- 188, 2001. 

Htay Aung Win. Model equation of surface waves of viscous fluid down an inclined plane. J. Math. Kyoto Univ., 
33(3):803-824, 1993. 

Jun Yu and Yi Yang. Evolution of small periodic disturbances into roll waves in channel flow with internal 
dissipation. Stud. Appl. Math., 111(1): 1—27, 2003. 

Kevin Zumbrun and Peter Howard. Pointwise semigroup methods and stability of viscous shock waves. Indiana 
Univ. Math. J., 47(3):741-871, 1998. 

Kevin Zumbrun. Planar stability criteria for viscous shock waves of systems with real viscosity. In Hyperbolic 
systems of balance laws, volume 1911 of Lecture Notes in Math., pages 229-326. Springer, Berlin, 2007. 

Kevin Zumbrun. Numerical error analysis for Evans function computations: a numerical gap lemma, centered- 
coordinate methods, and the unreasonable effectiveness of continuous orthogonalization. ArXiv e-prints, April 
2009. 


Brown University, Providence, RI 02912 
E-mail address'. Blake_Barker@brown.edu 

University of Kansas, Lawrence, KS 66045 
E-mail address', matjohn@ku.edu 

Institut de Mathematiques de Toulouse, Toulouse, France 
E-mail address: pascal.noble@math.univ-toulouse.fr 

Universite Lyon 1, Institut Camille Jordan, INRIA EP Kaliffe, Villeurbanne, France 
E-mail address: rodrigues@math.univ-lyonl.fr 

Indiana University, Bloomington, IN 47405 
E-mail address: kzumbrun@indiana.edu 


38 



